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ABSTRACT 

The Anderson model of wave localization in random media is invoked to study the effect of 
solar wind density turbulence on plasma processes associated with the solar type III radio burst. 
ISEE-3 satellite data indicate that a possible model for the type III process is the parametric decay 
of Langmuir waves excited by solar flare electron streams into daughter electromagnetic and ion- 
acoustic waves. The threshold for this instability, however, is much higher than observed Lang- 
muir wave levels because of rapid wave convection of the transverse electromagnetic daughter 
wave in the case where the solar wind is assumed homogeneous. Langmuir and transverse waves 
near critical density satisfy the loffe-Riegel criteria for wave localization in the solar wind with 
observed density fluctuations -1%. Numerical simulations of wave propagation in random media 
confirm the localization length predictions of Escande and Souillard for stationary density fluctua- 
tions. For mobile density fluctuations localized wave packets spread at the propagation velocity of 
the density fluctuations rather than the group velocity of the waves. 

Computer simulations using a linearized hybrid code show that an electron beam will excite 
localized Langmuir waves in a plasma with density turbulence. An action principle approach is 
used to develop a theory of non-linear wave processes when waves are localized, A theory of 
resonant particles diffusion by localized waves is developed to explain the saturation of the beam- 
plasma instability. It is argued that localization of electromagnetic waves will allow the instability 
threshold to be exceeded for the parametric decay discussed above. 
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1. Introduction 

In the usual picture of weak plasma turbulence theory, the elementary excitations are taken to 
be plane waves that propagate unimpeded to infinity. These elementary excitations are the standard 
plasma waves: Langmuir waves, ion-acoustic waves, electromagnetic waves, etc. {TsytovicK 1977), 
The plasma which supports the waves is assumed uniform, infinite and homogeneous except for the 
possible inclusion of a constant uniform magnetic field. Density gradients or fluctuations are 
treated as weak perturbations which refract waves smoothly or cause scattering of one elementary 
excitation into another. 

The picture that waves are freely propagating for times in the remote past and far future, and 
that interactions all take place in some intermediate time is a fundamental assumption in weak and 
multiple-scattering theories. On the other hand, a Langmuir or electromagnetic wave which ori- 
ginates and propagates near the critical density will be strongly scattered and modified by relatively 
small density fluctuations, and the scattering theory picture may not be valid. 

In his famous 1958 paper, P. W. Anderson investigated another picture for the elementary 
excitations in a randomly disordered system. In studying the absence of spin diffusion in silicon at 
low temperatures, he suggested that spin wave excitations, when the localized site frequencies are 
randomly distributed, are time-stationary wave packets localized in space with approximately 
exponential fall-off in amplitude away from the site of localization. The Hamiltonian he considered 
is equivalent to that of a wave excited in a medium with a random index of refraction, and so is 
equivalent to the problem of a wave in a plasma with a density that fluctuates randomly about some 
average density. 

In this thesis we adopt the point of view that the elementary excitations of a plasma with tur- 
bulent density variations are, in fact, weakly interacting localized (and therefore quantized) waves. 
The significance of this picture is profound. It suggests that a gentle bump-on-tail instability will 
excite localized modes with central wave numbers resonant with the beam. Therefore the amount 
of beam energy extracted should be greatly reduced in comparison with a beam propagating 
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through a uniform plasma. Wave-wave coupling processes will be modified. First, coupling con- 
stants must be reduced by a geometric factor that accounts for the spatial overlap between the cou- 
pled waves. Second, there is no wave convection for the localized waves, so the threshold for 3- 
wave processes is correspondingly modified. 

One significant difference between a turbulent plasma and a crystalline lattice with random 
impurities is that the former does not have a time stationary Hamiltonian. Our point of view is that 
the turbulence is externally or self-consistenily driven and independent of the waves in which we 
are interested. The effect of time dependent turbulence on localization is therefore an essential 
consideration. 

The motivation for this study is to attempt to clear up certain problems in the understanding 
of the type III solar radio burst. We adopt the model for this phenomenon presented by Cairns 
(1984) and Lin et al. (1986), which is, in brief, the following. High energy electrons, accelerated 
by solar flares, stream through the solar wind plasma and excite Langmuir waves through a gentle 
bump-on-tail instability. These Langmuir waves then decay parametrically into an lon-acousiic 
wave and an electromagnetic wave (at the local plasma frequency). These electromagnetic waves 
are then observed at the earth by microwave receivers as the narrowband emission that is charac- 
teristic of the type III burst. Although there is strong experimental evidence for this process, some 
theoretical issues remain unresolved. 

One remaining issue is the ability of the beam to propagate out to 1 AU without disruption 
while maintaining a well-defined positive slope. Another is that the electromagnetic decay is a 
convective instability with growth lengths larger than are consistent with the problem, A third is 
the spikey nature of the wave envelope of the Langmuir waves. It has been apparent for some time 
that the resolution of these difficulties must lie in density inhomogeneities in the solar wind plasma. 
While evidence for the existence of the density turbulence is strong, the details of the turbulence 
are still unclear. Efforts to understand the effect of the turbulence within the framework of tradi- 
tional weak turbulence theory have generally required very special assumptions about the nature of 
the density fluctuations and often more problems are raised than solved. Localization, on the other 
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hand, is a universal characteristic of random media (for sufficiently strong turbulence), as will be 
discussed, and any theory must consider the effects of this as well. 

The effect of density fluctuations on 3-wave processes has recently been discovered to have a 
pronounced effect in another context. Rose et aL (1987) have shown in plasma simulations that the 
ion-density fluctuations produced by Stimulated Brillouin Scattering, when sufficiently strong, will 
spatially localize the Langmuir waves produced by Stimulated Raman Scattering (SRS). Since the 
allowed Langmuir waves are now quantized, the SRS instability may or may not be suppressed, 
depending upon whether any of the quantized Langmuir waves satisfy the 3-wave frequency and 
wave-number matching criterion. In this case, however, the trapping of the Langmuir waves is 
within a single density depression and can be interpreted semi-classically, which will not be the 
case for the problem considered in this thesis. 

This thesis is organized as follows. In order to motivate this thesis and set the problem in 
context we begin in section 2 with a review of the type III solar radio burst problem, presenting 
essential satellite data and current models that attempt to explain these data. 

In section 3 we discuss the action principle in plasma theory which will provide the theoreti- 
cal framework for the remainder of the thesis. In section 4 wave propagation in random media is 
discussed. Multiple scattering theory is reviewed to explain why it is inadequate for this problem. 
The Anderson model of wave localization is then discussed. The results of 1-D numerical simula- 
tions of the time evolution of gaussian wave packets in random media are presented in section 5. 
These simulations verify the localization length predictions of Escande and Souillard (1984). The 
effect of mobile density fluctuations is also studied in secion 5. 

The bump-on-tail instability in random media is investigated in section 6. A theor>' for this 
instability is developed which takes into account that the excited waves are localized by the Ander- 
son mechanism. Results of a 1-D hybrid simulation are presented verifying this model. A theory 
for the saturation of the bump-on-tail instability by resonant diffusion of beam particles by local- 
ized waves is developed and compared with numerical simulation and experiment. 
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Section 7 begins with a discussion of 3-wave decay processes based on an action principle. 
Convective effects in homogeneous media and in media with a uniform density gradient are dis- 
cussed to conclude that the 3-wave process under consideration has too high a threshold. Previous 
work in 3-wave processes in inhomogeneous media is reviewed. Implications of Anderson locali- 
zation for the 3-wave process are then developed. 
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2. Type III Solar Radio Burst — Observations and Current Theory 

We begin with a discussion of the type III solar radio burst phenomenon and a presentation 
of current satellite data. The arguments in favor of a 3-wave decay model for the phenomenon are 
reviewed, and experimental measurements of density fluctuations in the solar wind are discussed. 

The type III solar radio burst is one of the oldest problems in radio astronomy; however, it is 
not yet understood in theoretical detail. Originally discovered by Wild and McCready (1950), the 
phenomenon is observed by an RF spectrograph at the earth as a narrowband RF emission, starting 
at -250 MHz and decreasing over a time period of -5-15 sec down to -5 MHz. Observations at 
lower frequencies using ground-based systems are not possible because of ionospheric reflection. 
Satellite observations have extended the lower range down to -10 kHz. Very often emission at or 
near the 2nd harmonic is observed simultaneously with the above. Wild (1950) proposed that the 
burst was caused by a disturbance which started low in the solar corona and propagated outward, 
exciting plasma oscillations of progressively decreasing frequency. Later it was established that the 
disturbance propagates at a velocity -c/3 from which it was inferred that the disturbance was 
caused by electron streams {Wild, 1954). More recent satellite observations, discussed below, have 
confirmed that the source of the burst is associated with electrons propagating along open magnetic 
field lines away from the sun. An interesting peculiarity is that there are some bursts whose fre- 
quency decreases and then increases again. These are interpreted as bursts due to electrons travel- 
ing on closed field lines associated with coronal loops and reluming to regions of higher plasma 
density near the sun. 

IMP-6 satellite observations at I AU {Fainberg and Stone, 1974) showed that the type III 
radio emission extends down to -30 kHz and the burst duration extends to many minutes or hours. 
Fainberg and Stone were also able to show that the emission region follows the Archimedean spiral 
of the open magnetic field lines predicted by Parker (1958). More recent satellite observations in 
connection with observations from ground-based solar observations have established that electrons 
with energies from a few keV to -100 keV, accelerated in solar flares, are responsible for type III 
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emission. 

The most detailed data on the type III solar burst come from the ISEE-3 satellite in orbit at 
the Lagrange point -258 Re upstream (towards the sun) from the earth. The satellite carries the 
University of California solar electron experiment, which provides electron energy measurements 
from 2 keV to 1 MeV, and the joint TRW-JPL-Iowa plasma wave instrument which provides elec- 
tric field measurements in frequency channels from 17.8 Hz to 100 kHz. Lin et al. (1981) were 
able for the first time to measure the electron distribution function parallel to the IMF (interplane- 
tary magnetic field) and show that it develops a well defined bump-on-tail distribution. This study 
was continued and results more fully interpreted by Lin et al. (1986), These results form the basis 
for the rest of this thesis and are discussed in detail here. For detailed history of t> pe III observa- 
tions with pertinent references see the review by Goldstein (1983). We shall next discuss the data 
from the Lin et al. (1986) study for a typical solar type III radio burst. 

Type ni Event — 11 March 1979 

Table 2.1 shows typical plasma parameters for the event. Figure 2.1 shows the high-energy 
non-thermal electron flux associated with the event along with electric field measurements for the 
II March 1979 event. The emission in the 100 kHz channel starting at -10:40 UT is characteristic 
of the type III emission. It is interpreted as radio emission upstream of the spacecraft which pro- 
pagates to the spacecraft. Onset of emission at lower frequencies occurs later. Activity in the 17.8 
kHz channel, corresponding to the local plasma frequency, occurs with the arrival of electrons in 
the 8.5 keV channel and corresponds to times when a positive slope bump-on-tail is observed in the 
reduced parallel electron distribution function shown in Figure 2.2. 

Figure 2.3 shows a portion of the same event in high time resolution. The important points 
to observe are the intense spikiness of the Langmuir waves in the 17.8 kHz channel, and the coin- 
cidence of activity in the 100 Hz channel with the most intense spikes of the Langmuir waves. 
Early on there was some ambiguity about whether the radiation observed by satellite is fundamental 


or 2nd harmonic emission, but it seems certain from these data that the emission is fundamental 
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emission. 

The foUowing assumptions are made in interpreting the above results. The dominant wave 
number for the Langmuir waves is given approximately by the beam resonant wave number 
where is the beam velocity. The ion-acoustic waves observed in the 100 Hz channel 
are actually much lower in frequency and are Doppler shifted by convection past the satellite by the 
solar wind, i.e., <o' = V^^f/COs9 (0 is the angle between the solar wind direction and the inter- 
planetary magnetic held). Assuming one calculates the wave number for the ion-acoustic 

waves as 


, ^observed 

K = 

V COSd 

The actual frequency is then inferred from 

Based on the above data Cairns (1984) was the first to propose that the important emission 
mechanism for the type III emission is the parametric decay process L (Langmuir) — > T (transverse 
electromagnetic) + S (ion-acoustic). The argument was based on considering the volume emissivity 
(brightness temperature) of various 3-wave processes. The weak turbulence process L — ^ 7 + S 
(random phase approximation) produces emissiviiies that are far loo small. Further evidence for 
this process was observed in the experiment of Whelan and Stenzel (1985). In this experiment 
Langmuir waves were excited with a weak beam. Electromagnetic radiation at the plasma fre- 
quency was observed to occur with growth that closely followed the growth of ion-acoustic waves. 
The ion-acoustic waves were observed to have a wave number k^ = ki, the wave number of the 
Langmuir waves. This process was also observed in a 2-D numerical simulation {Pritchett and 
Dawson, 1976); however, they interpreted it as a two-step process: 

L -^L' + 5 
followed by 

In space, however, the threshold for the process L^L' -\-S is not satisfied by several orders of 
magnitude {Lin et al, 1986). 
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A fundamental difficulty with this interpretation is that in homogeneous media the process is 
in fact a convective instability, since the threshold for absolute instability is not exceeded as a 
result of rapid convection of transverse wave energy and the high damping rate of the lon-acoustic 
waves. This is an involved argument and will be taken up again in section 7.2. In fact, the solar 
wind is not uniform and we will argue that the solar wind density turbulence is fundamental for 
understanding the type III process. 

We note briefly that the frequency and wave matching coupling conditions will give a 
daughter ion-acoustic wave with wave number = The transverse wave will be polan/ed with 
the electric field along the direction of the Langmuir pump and propagate transverse to this direc- 
tion with wave number 


(Lin et aL, 1986). We note that 

kj ^ and + (o^ = . 

This is discussed fully in section 7. For the March 1 1 event this gives A 7 = 360 km. 

Many efforts have been made to understand the type III mechanism as the result of a modula- 
tional instability leading to solilon formation and collapse. The experimental evidence does not 
support this, and furthermore the threshold condition 
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^ 0 - 




innkT 

where AA* is the half width at half maximum of the spectrum of beam excited Langmuir waves is 
not satisfied. This is a complex subject with a very extensive literature, and shall not be discussed 


further here (but see Lin et aL, 1986 and references therein). 
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Solar Wind Plasma Density Fluctuations 

Figure 2.4 shows a 3-hour time averaged power spectrum of solar wind density fluctuations 
that was obtained by Celnikier et at, (1981) by measuring the phase shift in a radio beam propagat- 
ing between the two satellites ISEE-1 and -2. While fluctuations occur at all wavelengths, it is par- 
ticularly noteworthy that fluctuations with wavelengths of -100 km have 6n/n>\%. Little is 
known in detail about these fluctuations, and there are no results available simultaneous with a type 
III burst. On the other hand, density fluctuations of this level appear to be a universal feature of 
the solar wind, confirmed both by scintillation measurements {Coles and Harmon. 1978) and by 
recent Giotto satellite data (/?. P. Lin, unpublished data). The fluctuations seem to be nearly isotro- 
pic with a ratio of fluctuations parallel to the interplanetary magnetic field to those perpendicular 
to the magnetic field being less than 2:1. Although these fluctuations may be a kinetic drift wave 
instability {Mikhailovsky, 1983) resulting from initial solar wind inhomogeneiiies at the corona, 
magnetic field fluctuations suggest that they might be electromagnetic in origin {M. Goldman, 
private communication). 

For some time now there has been an awareness that density fluctuations must play a role in 
the type III problem. The first problem that has bothered people is that Langmuir waves come in 
clumps or packets -70 km in extent (see Figure 2.3). Since the electromagnetic radiation is very 
uniform in time (i.e.. Figure 2.1, 100-31.6 kHz channels) one concludes that this dumpiness is spa- 
tial rather than temporal. Smith and Simes (1977) studied the propagation of Langmuir waves 
through a turbulent plasma using ray-tracing techniques to conclude that enhancements and reduc- 
tions would result due to random ray focusing. They concluded, however, that the main cause of 
the clumping was that the rays would be resonant with the beam leading to growth in certain 
regions. Implications of this reasoning for 3-wave processes have been examined by Melrose et al. 
(1986). The remainder of this thesis is an effort to understand the effect of density turbulence in 
detail. In particular we will argue that the correct effects of the density turbulence can only be 
predicted by studying the full wave equation in random media, and that ray tracing is not a satis- 
factory approximation when attempting to understand wave behavior near critical density. 
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It should be noted that much of the effort of solar radio astronomers has employed ray tracing 
techniques to map observations at the earth back to the source region in the solar wind. Although 
much of this work is probably valid when identifying source regions, conclusions about the exact 
size or geometry of source regions is questionable. 
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3, Theoretical Introduction 

In order to make theoretical progress with this problem we use an action principle approach 
for a wave-plasma system {Dewar^ 1970, 1972,1973; Johnston and Kaufman^ 1978a,Z>, 1979; Kauf- 
man et al., 1984; Crawford et al.y 1986; Similon et aly 1986; Kaufmany 1982, 1987). In this 
approach we begin with a many particle Lagrangian 

l = (3.1) 

"2c 8;r 

where ^ represents a sum over all particles in the system. 

i 

This may be re-written, using a canonical transformation, as 

^ = X [p-^i - (-V,. a; At,))] - j [l £ I - - 1 fi I -] (3.2) 

where 

( 

I ^1 

/i, = zr- Pi A (,r, ) +ej(t> (x, ) 

2 m c 

\ / 

is the single particle Hamiltonian. The least action principle with this IV -body Lagrangian, in prin- 
ciple, determines the time evolution of any (non-reiativistic) system of interacting particles and 
electromagnetic fields. 

We shall divide the particles into groups according to the following distribution functions: 

1) /o^(.v,v,r) = non-uniform solar wind background electron distribution; 

-) /oi('Y,v,0 = non-uniform solar wind background ion distribution; 

3) /i,e('V,v,f) = electron beam distribution function. 

In the particular type III problem, the beam density is 5-6 orders of magnitude below the 
solar wind density, and therefore, we make the following separation of terms. We define 

and to be self-consistent fields arising from strong turbulence associated with non- 
uniformities in the solar wind, and then define = E - and 
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= 0 — (p^^. 


We use these definitions to expand the Lagrangian as: 


^=I 


<BO 


2 2 

Pi -x. - -^+-^Pi T-S- IA^‘^(jTi)l-ei P^^{X.) 

Irrii nXiC ImiC' 


Zni 


iBO 




m,c 


IrriiC^ 


- [l£<‘>|-- IB<')|-j 


+ y-^ aC)(;c..)-a5C(.v,.) 

i rtijC’ on-’’ 


(3.3) 


O' l.^"'(.V, )! = -., »i, 

/W/T 2w,6’^ 


<6 


_fi 


/WyC' 


where /*bg refers to background particles and refers to beam particles, and / =/bg + 

Let us begin to investigate equation 3.3. We choose our definition of the self-consistent 
fields such that in the absence of beam particles, the first line completely describes the evolution of 
the solar wind background particles and their self-consistent fields, i.e., variations of the first line 
with respect to 0^^, jtbg. and pbg vanish. In general, for density distributions and magnetic 
fields present in the solar wind, this is an intractable and unsolved problem in strong turbulence 
theory which will produce a whole zoo of drift waves, magneto-acoustic waves, etc. We do not 
discuss this pan further, but imagine the problem solved, giving us the turbulently fluctuating back- 
ground distribution functions and self-consistent fields /oiU, v,r), /o^(.r,v,r), A^^Cr,/), and 
(p^^(Xyi). We shall be concerned with the further evolution of the system due to other waves and 
the beam particles, having assumed that we can separate out the strong turbulence portion of the 
problem (not amenable to an order by order perturbation treatment). 
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In typical approaches to such problems the Lagrangiun is divided into a time averaged part 
and a fluctuating part. The time averaged part is defined as 



and the fluctuating part as 


L -L . 

A suitable canonical transformation will transform the fluctuating part to higher order in the expan- 
sion, as discussed below. The order of a term in the expansion is given by the sum of the powers 
of the field quantities. In the approach developed by Johnston and Kaufman (1978), this transfor- 
mation is accomplished by a Lie transform on the particle coordinates r,, p,, giving a time averaged 
(guiding center) part plus a higher order term that will involve the electromagnetic fields. The 
Lagrangian, expressed in the new variables, is then solved, order by order, using a least action prin- 
ciple. 

Thus we define the transform of a function on phase space (.v,p) as (Cary, 1978; Lichtenberg 
and Lieberman, 1983): 


e^ = \ + L + \L-+jL^+... 
where, it A = A (.x,p ), then 


i 


dw 3A 
3AT, dp, 


dw 3A 

dpi d.r, 


s { w,A 1 


for any given w = w(r,/>). 

In this transformation scheme the new coordinate variables become 

P=e‘-p 
X = e'-.x 

and the particle Hamiltonian is 
K = e‘-H 
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if H is ordered, 

then we define terms in the transformed Hamiltonian: 


ATW = //<"> 

^(I) = //0) + £^(0) 

A:(2)=//a)+£i/(i)+l£2//(0) (3 4) 

AT<^)= + LH^-\ etc. 

L 2 

In our problem we shall not consider terms over 3rd order. We then have the action 




-A:, 


V J 

Where is independent of the particle coordinates and therefore is unchanged by the Lie 

transform. This transforms to 


I 

We choose w such that vanishes, or = vi^//^}. Since, in fact, { w, //^} = for 

dt 

unperturbed trajectories, w = -j integrated over unperturbed particle trajectories 

x = — t; P = Pq^ 

m 

The third line in equation 3.3 will time average to zero (since we assume that the turbulent 
spectrum means that there is little self-consistent field energy at a frequency and wave number 
resonant with the wave under consideration. This is supported by I. Cairns s (1984) argument that 
random phase processes do not explain the intensities observed. We shall also assume that the last 
line in equation 3.3 can be ignored. This is equivalent to assuming that the turbulent fields have 
little energy resonant with the beam electrons, and so therefore the effect on beam propagation can 
be neglected. When the Lie transform is applied we will get a term of the form 
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where F is a function linear in its arguments. This term is very small compared with those of the 
2nd line also quadratic in and therefore can be ignored. We shall from here on ignore ions in 

the quantities involving wave-coupling terms since the coupling terms go as — and — » — . 

tn nig trii 

Wave Propagation in a Turbulent Plasma 

With the beam switched off, the evolution of a wave of sufficiently small amplitude that its 
effect is negligible on the development of the strong density turbulence, will be governed by the 
second line of equation 3.3: 

5^’>= f'' dt Pi At'Hxi) - -At I I - - a I ' ' - ' ^ ) 1 

•'^0 ■ niiC Om.n- 8;r ^ I 

I iflo * 


IrriiC 


If we take 


— Pi • aP (.y,) - e, 
rriiC 


//(■)=£ 
i 

then w = -J dt {Johnston and Kaufman, 1979). 
Using this, we find 




and 


S = f '■ dt I { w, - £ -At I aC) (.y.) 1 2 - ^ J d^Y 1 £(') 1 2 - 1 bO) 1 2 

Every term in this expression is quadratic in the fields. If we can use the eikonal representa- 
tion of a wave, then 

A = A + 

where 0 is a rapidly varying phase and A(x,t) an amplitude that is slowly varying in time and 
space. McDonald et aL (1985) have shown that this is equivalent to 


S = j dt co \ x,t) ' D {k,o>\x,t) 'A{k,co\ x,t) 
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where D is the normal dispersion tensor. Variation with respect to the held amplitudes then gives 


5S 

3A 


= 0= D (A’,co) A(jt,a)), 


which is the usual dispersion relation. 

For example, if we only consider electrostatic waves of the form 


= 2 e'(^-“') + cc. 


then 




i(kXj-mt) 


kVj — CO 


+ CC 


and 


S=jdtj 

= _ f y 1 iL ^ I . 

J '“2 m; (^-v, - co)- 8;r J 

Using the Vlasov representation of panicle phase space density: 

X = Z J = I V, t) 

i i 

we find 




The quantity in brackets will be recognized as e{k,ay), familiar from the usual plasma kinetic 

I 

I 

theory (Lifshitz and Pitaevskii, 1981). It should be recalled in examining equation 3.4 that this 

i 

does involve the turbulent particle motions, since /(.v, v, t) involves the instantaneous panicle posi- 


tions and velocities. 


While the Lagrangian approach does lead nicely to a dispersion relation in A: -space in the uni- 


form density case, it does not immediately give us a convenient equation for the time dependence 


of wave amplitude in the turbulent case. Therefore, we back up a step and redo this part of the 


problem using fluid equations of motion. {Nicholson, p. 27ff). Actually, the fluid equations of 


motion follow from a Lagrangian such as equation 3.2 by taking suitable moments of the 


d^y 


4ne“ f (.V, \\ i ) 
m (jb’-o>r 


- 1 


(3.4) 


I 
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distribution function and expressing them as fluid quantities (Dewar, 1970). For the case of longi- 
tudinal waves: 


9/7- 

IT"-'' ■(”•'■) 

St m mn^ 

V E =-4;re(/7^-/7/) 


We will further assume that rii = /7/(.v, v, t) is already known from the strong turbulence prob- 
lem and that quasi-neutrality holds. This latter condition is quite good in the type III problem, 
since kp - 22 m, while ^Langmuir 
We get 


dr- 


(V ‘£) = V 


^ ^ ykT ^-E 
rn dx- 


m 


“4;r e V *([V -(n^v)] v) . 


The last term is O(v') and will therefore be ignored. Therefore, for longitudinal waves 


^ +3vi-j^ 

where = 3v,^. 
m 

Beginning with Maxwell’s equations, it is elementary to derive a similar equation for elec 
tromagnetic waves 


9“£ ^ d~E 

') 47Te* 

where again co“(x)s n,(.\:,r) is determined by the strong turbulence problem. 

^ m 

The properties of these wave equations form the content of the next chapter in Anderson 
localization. Upon concluding with that, we will return in chapters 6 and 7 to our Lagrangian to 
understand the beam-plasma and wave-wave interactions. 


• We discuss, finally, our assumptions about the geometry of this problem and the efrect of 
magnetic fields. At 1 AU, the IMF has an average value of 5y (ly = I x 10”^ gauss). An electron 
plasma wave perpendicular to the magnetic field will have a frequency at the upper hybrid 


18 


frequency. 

= 1 .0001 16 CO' . 

As this is much smaller than the fluctuation in the local plasma frequency (-1-2^), due to density 
fluctuations, we will ignore this effect. The major effect of the magnetic field will be to constrain 
the (weak) electron beam to flow along the magnetic field lines, and therefore the wave vector of 
the Langmuir waves is directed approximately along the magnetic field. 
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4. Wave Propagation in Random Media 

Introduction 

A wave with wavelength A, propagating in a non-dissipadve three-dimensional medium with 
a low density of scatterers with a scattering mean free path I* such that will exhibit 

diffusive behavior when observed at a distance L':^ I*. As the density and strength of the scatterers 
is increased the diffusivity will decrease until, at some value of diffusion will cease alto- 

gether according to the theory of Anderson localization. The wave will then be localized in space 
in a time-stationary quantized eigenstate. 

Wave propagation in a collisionless turbulent plasma such as the solar wind may provide one 
of the few observable examples of strong wave localization in a classical system. In general it is 
difficult to find systems which satisfy the loffe-Riegel criterion for localization, namely that the 
scattering mean free path should be of the same order of magnitude as the wavelength and that the 
dissipation length be much longer than the mean free path /* {Anderson, 1985). In a weakly coUi- 
sional plasma, such as the solar wind, the dissipation rate y for the electromagnetic waves will be 

of the order of — ^ - 10“'^ and for the Langmuir waves of interest (since A » Ar,) — - 10"^^. 

o);, 

Another factor that makes plasma waves attractive for studying localization is that the dielectric 
Afi 

contrast due to density fluctuations, where e is the dielectric constant, can be very large, since 
for plasma waves the ensemble average <e> s 0 and for electromagnetic waves near the plasma fre- 

quency e — is very small. 

The condition A - / * is in fact a weak condition, and values of I* IX as large as 50 have been 
observed to produce localization as argued by Trawal et al, (1986), who attributed the observed 
vanishing of positron mobility in gaseous helium to localization effects. The ratio /*/A appears to 
be within the appropriate range of values for Langmuir and electromagnetic waves associated with 


the type in solar radio burst. This is in contrast to the failure to observe strong localization effects 
in the propagation of a laser beam through a slurry of small dielectric spheres {Watson et al„ 
1986). The problem in this experiment was that the dielectric contrast was too low for latex 
spheres, and so the mean free path was too large. 

Before discussing Anderson localization, we begin first with a review of multiple-scattering 
theory in order to develop some of the language which will be important in localization theory and 
to show how multiple scattering theory breaks down in the regime where localization is important. 


4.1. Multiple Scattering Theory 

In discussing multiple scattering theory we rely on standard Green’s function techniques dis- 
cussed, for example, in Abrikosov et al. (1975) and applied specifically to the problem of propaga- 
tion in random media in Lax (1981) and Rickayzen (1980). 

Multiple scattering theory begins with a wave equation of the form 

= V(x) v(x), (4.1.1) 

One imagines that as r -» ±<» the wave is in a packet which, however extensive, in a region where 
K(jr)-^0, i.e., no scattering takes place. In this region the wave will asymptotically, as t 
be of a form 






where a>j; = a^ + k^c~. One uses perturbation theory to develop solutions which will again be of 
the same form for r -» «>. A comparison between y(t -®®), and y{t -> oo) gives the scattering 
amplitude. 

If we Fourier transform equation 4.1.1, we get 


[o>2 - v(k ) = Jv(k-k')v(k')dk'. 


Let us define a Green’s function G{x,x', t,t') such that 
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The significance of ihe Green’s function is that if we have some initial condition, 
VTofjr, / = 0) , then 

VU.Os Jd’.vG(jr,ar',/,0) VToU.' =0) 

is a solution of the wave equation. We also define Go{x,x\ t, t') such that 

G^{x-x',t-t') = S{x-x')S{t-n 

In Fourier space this becomes 


a- , a- 


hx'^ 


[ti}2 - - c'k^] Goik, co) = 1 

or 


Goik, CO) = i . 

'y T I > 

0>“ - 0>“ - C'k * 

It can be shown that 

G (k, (o) = Go (k,<o) + jdk' Go(k,co) v (k -k')G (co, k ') . 
This can be iterated to for a perturbation expansion in v to obtain 


G = Go + GovGo + GovGo\Go + • . . 


where Go and v are integral operators (or matrices). 

One represents G diagramatically as => and Gq diagramatically as — >. 

We can represent 4.1.2 as 


u ^ J 


h. 1^. 




H.-VL k 






w w • k 


where x represents the interaction v. One defines a “self energy” ^ {k, co) with graph 


(4.1.2) 



) 


such that 


G (k ,o>) = Go(k,to) + Go (k ,co) ^ (A- ,a>) G„ (A ,o>) + Go(k) ^(k) Go(k ) + . . . 
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This can be solved as 


G(k,co) = 


Go'' 


In order to see how this theory works we consider a simple one-dimensional example. Let us 
consider a single perturbation in an otherwise uniform medium of the form 


Vo -lJ 2 ^x<lJ 2 
Q otherwise 

the diagrams for the self-energy are 








+ - 




plus terms of higher order in v. 

Using the rules in the above cited references, the first order diagram for the self-energy just 
leads to a shift in the average plasma frequency. The second order diagram in the self-energy may 
be evaluated as: 


1 v(^)l-Go(A -<?,co) 


where - — corresponds to the normalization of the inverse Fourier transform, and where 
27 tL 


2v’o 

’ \ ~ J 


Therefore, 


. 1 (•> sin-(r//^/ 2 ) 

\(.k , (o) = ^~ ] d(/ — p 

) CO- - oy- - c-{k - qY 


and then 


Gik^co)-- 


y~ - co^- c~k~ - ^ (f’, co) 


The pole avoidance procedure is that a> -> co -i- /t] corresponding to the retarded propagator. For Vo 
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sufficiently small, we will only be interested in ^ when o)^-<oj-c^k^ = 0. In this case with the 


2/2 


approximation ^(Jk, to)= ^ we find 


G(k,a>)s 


„2_,.,2_^2t2. * ^0 // 


- 


2kLc^ 

And we see that the self energy then represents a lifetime for the excitation. 


For a large system with many scatterers we replace — with N, the volume density of scatter- 
ers, and we take an ensemble average. Then 


<G (k, a>)>s 


1 


2 2 2 , 2 i^^olc 

(a c k 


2c-k 

In order to ensure that higher order terms in the expansion may be neglected in the self-energy, we 
require that But by examining the third diagram in '^{k, <o) one estimates that 


'( 3 ) 




’( 2 ) 


Vo V 

Therefore, we expect that a perturbation theory will break down when — ; 1. As we shall see 

in the next section with vo=5o>^, this is the regime when Anderson localization becomes important, 
and will be the parameter regime of importance for the type in problem. 

If we have an initial amplitude i/(ato, /o)» we can compute the transport of probability density 


as 


Using the Green’s function language 

I yf(x, f) 1^ = v(x,t) v /’( at ,/) = I d^Xo d^x’oG(x, t,Xo, to) Vixo) V (x'o)G\x, t,x'o, to) (4.1.3a) 
The diagram for this is 
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(4.1.3b) 


where = G= s=> 

is the retarded Green’s function and 


is the advanced Green’s function. 

Expanding (4.1.3) gives diagrams which include the "ladder” diagrams: 


X 

I 




and the maximally crossed diagrams: 


*>1 X 


X , A 
t 


Langer and Neal (1966) showed that in evaluating diffusivity in random media these two sets of 
diagrams axe the most important^ and thau in fact, their sum diverges unless one adds in finite life- 
times in the unperturbed propagators (the Go or light lines in the diagram). 


Particularly important is that the ladder and crossed diagrams interfere constructively for 
exact backscattering. This is a consequence of the time reversal symmetry of the backscanering 
processes (Gor*kov, 1977). Because this constructive interference is a symmetry that survives 
ensemble averaging of the above two-particle diagrams it is expected to be a dominant effect and 
lead to enhanced backscattering. This enhancement has been observed recently by Van Albeda and 
Lagendijk (1986) in an experiment that measured the backscatter of laser light from a slurry of 
glass beads, and is also important in understanding experiments on electron mobility in thin two- 
dimensional films. The parameter regime where enhanced backscatter is observable, but behavior is 
otherwise diffusive, is often referred to as “weak- localization.” 
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One of the major difficulties in studying localization effects analytically is that 
<G'* G^>5*<G '*><(?*>. Correlation effects are the essence of Anderson localization theory. 
Enhanced backscatter, for example, depends upon these correlations. We shall return to the issue 
of two particle Green’s functions later in the next section. 

Before concluding this section on multiple scattering theories, we mention the effort of 
Muschietti et al. (1985) to understand the effect of 3-D density fluctuations on beam-plasma insta- 
bilities by considering the Fokker-Plank equation 

d, w* = 2% w* + dtj Dj, dij Wi + 5* 

with 

Dj, = c{q)qj q, 6(2k ■ q ) , 

where y* is the Langmuir growth rate for wave vector k, is the wave energy density, c{q) is the 
spectral density of fluctuations of wave number and is the rate of spontaneous emission. It 
was concluded that the diffusion of Langmuir waves led to the quenching of the beam-plasma pro- 
cess in the case of electron beams in the solar wind unless special assumptions about the anisotropy 
of solar wind density fluctuations are made. This theory is equivalent to a first order scanering 
theory, and ignores the possibility that as a result of infinite order scattering effects, a wave may 
remain in resonance with the beam for a long time. 

4.2. Anderson Localization 

As mentioned in the introduction, for sufficiently strong random scatterers waves do not 
diffuse in space, but rather are localized into normalized eigenmodes. The theory of wave localiza- 
tion in random media was first proposed by P. W, Anderson (1958) to explain the experimental 
result that electron spin polarization diffuses very slowly in highly doped silicon at low tempera- 
tures. Anderson studied the model Hamiltonian 

iaj = Ej + X (4.2.1) 


where a,(r) represents the quantum mechanical amplitude at the site i on a lattice, is a (weak) 
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coupling between sites and is the local site energy. The crucial point in the Anderson model is 
that Ej is taken to be a random variable, uncorrelated from site to site. For most of our discussion 
the probability distribution for £ is: 


t 


p{E)^ 


W 

0 


-WIl^E^WIl 

otherwise 


V 

although it is argued that the exact nature of the distribution is of little consequence, the width 
being the important factor. 


Equation (4.2.1) can be interpreted as the Hamiltonian for an election in a random lattice in 
the tight binding model, and therefore has been much studied in the context of solid state theory in 
order, for example, to understand the low conductivity of amorphous materials. It is relevant to our 
discussion since it is equivalent to the wave equation for either Langmuir or electromagnetic waves 
expressed as a finite difference equation. 


Localization in One Dimension 

Let us start with the wave equation 

-j^ = o>^(x)rp-v^V^-Y 

where 


for electromagnetic waves 



3v^ for Langmuir waves 


and look for time stationary solutions - 
in X where I is some small length: 





We rewrite this as a finite difference equation 


a>^\y = CDf(x) iVix + 1) + xpix - 1) - 2y(jr )) 

With the following definitions; 


(4.2.2) 


Vj = V(jl) 
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Vj = VO+i - Vj 

and ScOp. =a>p(xj)-<^^ where is the square of the average plasma frequency, equation 4.2.2 
becomes 

(o^Vj = ^ ) VO - -J7 (VO+I - VO) (4.2.3) 

In order to understand the nature of localization in one-dimension, we reproduce an argu- 
ment, originally due to Mott and Twose (1961) and recast by Escande and Souillard (1986). 

The general rule for solving a second order finite difference equation of the form 4.2.3 is to 
use a leap-frog method, by which i/' and y/ are alternately advanced. This gives a solution that is 
accurate to 0{^) rather than just 0(Ar) (Birdsall and Langdon, 1985). 

This gives first: 

-£^(1 - -^)Vy = - V) (4.2.4) 

and second: 

VO+i * V0>i - Vj (4.2.5) 

where we define 

V 

Vj 

E ' k-v- 

and is strictly a definition and its interpretation as a wave number will be dis- 

cussed later. 

We write (4.2.4) as 


Vj 


' 1 o' 

/ 

Vj 

vo>i , 

j 


1 

' / 



and (4.2.5) as 
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/ ^ 
vo>i 


f \ 

1 1 



\ y 


0 1 


V'/'+i 

/ 


which gives 




•i-E, |1 


f \ 

Vj 


/ 

¥j 

VO+1 

> 


, ‘J 


Vj 
\ > 

^Mj 

Vj 

\ / 


In the case where £^ = £, a constant, for all 7 , Mj is easily seen to be a unimodular matrix with 
eigenvalues 1 ± / 0(^)* Iteration in this case gives: 


A/,., ...A/o=A/" 

with eigenvalues , and therefore reproduces the solution expected for a uniform 

medium* 


In the case where the £/s and therefore the A//s are random an initial condition. 


( \ 
¥0 

W) 

V / 


, IS 


propagated as 


f \ 


V J 


= Afy A/y_2 ... A/i A/o 
We can define a Liapunov exponent 


f \ 
V'o 

Vo 
V ; 


Xs lim ^Tr(MMMj...M^Mo) 

It was first proposed by Mott and Twose (1961) and later proven rigorously by H, Furstenberg 
(1963) that for such a set of random matrices A >0 almost everywhere, i.e., except for £ taking 
values on a countable set {£,}. The significance of this is that almost any initial condition will 

developed an exponentially growing wave envelope going as at site 7 . An exponentially 

decreasing sequence is developed as follows. If 

II V^yll lIVyoH 

then we define 


1 


\ 



Vj 

i>; 
V y 
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If this is propagated backwards, 


V,.' 


vOo 

\ j 


/■a ^ 


*Af-' A/-> ...Af-' 




will be an exponentially decreasing sequence, although if continued beyond would again 
increase exponentially except for certain discrete values of E, We produce exponentially localized 
states by matching solutions from the left and right following an argument due to Mott and Twose 


(1961). Begin with two initial conditions, say. 


and 


( \ 
V^o 




and 


( \ 

1 1 

( \ 

Vn 

11 

Vo 

Vn 
V ) 


V'o 
^ / 


\ J 


, where n is large. One then finds 


( \ 


\ } 






one then adjusts E and 

Mfln 

/ 

such that 



f \ 

Vn 


J 




Vn 

\ / 


. Heuristically, one expects that this pro- 

r-i r-j 

duces a discrete set of values for E . 

The theory of localized states and the nature of the spectrum in one dimension was put on a 
rigorous mathematical basis by Kunz and Souillard (1980). An important point from localization 
theory is that in one dimension all states are localized. There are no good analytic expressions for 
localization lengths for WIV >\. Escande and Souillard (1985) determined the Liapunov 
exponents numerically for the product of random matrices The results are 

shown in Figure 4.1. The relevant parameters are 


WIE = SmpkWh 


and 
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The Liapunov exponent is 

Ic 

r=T 

where 4 is the localization length, and the “correlation length” is the distance over which o>^(.r) 
is constant. 

We now consider the correct interpretation of the quantity 
= (o>’ - (opfv . 

Clearly, when Vj^O for all j, as we have shown, k is the wave number of a plane wave propagat- 
ing in the medium. In the case where /:»1/^, i.e„ where there are many wavelengths (for the 
homogeneous case) in a localization length we would expect that the wavelength (measured as 
twice the average spacing between points where v^ = 0, say) would remain the same, since the 
fluctuating value of the spacing should average to zero, as should be the case when WIE is small. 
This is not entirely true, since although the fluctuations are uncorrelaied, the localized waves are 
correlated with the fluctuations, but this interpretation is borne out by the simulations in section 5. 

The physics of localization in a plasma is as follows for a plasma wave, or for an electromag- 
netic wave traveling near the critical density. The plasma may be viewed as a collection of weakly 

coupled harmonic oscillators. Each oscillator has the resonant frequency co^(.v) = + 5 co^(.r). If 

an oscillator is initially set in motion, it will couple to adjacent oscillators. Physically, how well 
this energy is coupled to adjacent oscillators will depend on two factors: 

(1) The dephasing rate adjacent oscillators: 

(2) The coupling strength -k-v~ which is a which measure of how well one oscillator can 
overcome the dephasing tendency of an adjacent oscillator, 

6(oz 

and thus the importance of the parameter ^ . 

k’v~ 

Thc other important parameter in localization theory is 
£ 
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which measures how near the frequency is to the “band edge/’ This will be discussed further 
below. 

Another way of viewing localization is that it is an infinite order scattering process in which 
a localized wave is self-consistently scattered back onto itself while the phases of the waves scat- 
tered in the “forward” direction are random and cancel out when viewed far in the forward direc- 
tion. 

Kunz and Souillard (1980) have shown that the spectrum of the equation 

V/ + [VA,>| + - 2 ^,] 

is 

[0, 4] + supp r 

where 

suppr^ { V'o IProb[V'£:{l/o“^t V"o + e}]>0 for all e>0} 

This is equivalent to the set 

{ o>2 I = cof + Scof] + [0, 4 1-//“] } 

for equation 4.2.3. For waves in a continuous medium, /->0 and the upper limit of this set has no 
meaning. 

This tells us nothing about the probability of finding a state located within a finite region of 
size much greater than the localization length. For example, for weak perturbations, a wave with 

approximate wave number k will have an energy which is near the value = cor - -l- 

with very small probability, since the probability of almost every wave peak aligning itself with a 
density minimum is very small. 

In localization theory an important parameter is how close co” is to the “band edge,” 

The concept “band center” does not have strict meaning since there is no upper band limit as is 
the case for a crystalline lattice. By making sufficiently large, we can make the localization 
length in one dimension arbitrarily large provided the condition remains met for plasma 


waves. 
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Localization in Higher Dimensions 

So far we have only considered localization effects in one dimension, and higher dimen- 
sionality effects must be considered. Anderson’s original approach was to take the Laplace 
transform of (4.2.1) which was written as 


|5( 


where the initial condition is ay(/=0)- Solving for /q gives 




'0* 


1 77 ^ 


/s-£* f /s-£* is-Ef 

\ / 

The question whether energy remains localized is equivalent to the question whether lim oqU) > e 

t—¥0O 

for some £>0. This is equivalent to lim fo{s)>e. This will in fact be the case provided that 
the self energy at site i = 0 


5^(0)=y + X — 1^5iZ«Zi2_ + , 

y is-£* ^ (/s-£*)(is-£,) 


(4.2.6) 

converges. Anderson studied this by considering the probability of the value of very long 
sequences when the £/s are random variables. Anderson argued that the self-energy would con- 
verge with probability 1 whenever WIV>{WIV)c, where {W/V)^ is to be determined. Anderson’s 
estimate was that (WI2V)c-2K ]n(WI2V)c where K is the connectivity of the lattice (from perco- 
lation theory) and is approximately z - 1.5 where z is the number of nearest neighbors. Thouless 
(1979) argued that this estimate is high because Anderson assumed that the contribution of each 
sequence was independent. Accordingly, a better estimate is 


K 

V 


= 4K 


A more transparent, but less rigorous argument for this is the following: The amplitude on a 
neighboring site will be 


aj~ 


V 
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Since the values in the denominator can range between ± W!2, a typical value is lV/4. If there are 
r neighbors then we require 


V \ W 

— = <— > or — >4Z. 


W/4 Z V 

Even this estimate, it turns out, is high. The most convincing numerical work on this was 
done by Vos/ii/to and Okazaki (1977), who solved the Hamiltonian Ea^ =£*, X where the 


j 

sum is over nearest neighbors, by diagonalizing the Hamiltonian on a 100 x 100 grid. They found 
that for W IV >6.5 a transition occurred between extended and strongly localized states. Thouless 
(1979) later argued that when W IV- 6.5 the state is localized with a localization length on the order 
of the size of the system. 


Actually, it is believed that VV/V=6.5 is a transition point between strongly and weakly 
localized waves in 2-D. According to scaling arguments (Abraiuims et aL, 1979) all states in 1-D 
and 2-D are localized, and for 3-D there is a critical value, {WIV)^ above which states are local- 
ized. Licciardello and Thouless (1978), looking at the sensitivity to changes in boundary values, 
concluded that 3-D wave functions are apparently localized for the diamond lattice for 
9<(W/V)^< 12 but the localization is slower than exponential. For the cubic lattice, localization 
occurs for W!V = 15. This was confirmed in a numerical study by McKinnon (1985) who studied 
the scaling with M of the localization length of long 1 -dimensional chains when adjacent 
chains were coupled together using nearest-neighbor coupling. His localization estimates were 
again based on a Liapunov exponent estimate. Also of great importance is that the normalized 
localization length h^^{M)IM in fact was found to follow the behavior \!M for strong locali- 
zation {WIV >25), thus suggesting that for strong localization, the localization length is only 
weakly dependent on dimensionality. 

A significant difference between our wave equation and that of the Anderson theory is that in 
our case there is no natural discrete lattice. In a crystalline lattice there is an upper limit on the 
wave vector in the unperturbed case corresponding to the first Brillouin zone. For Langmuir waves 
or electromagnetic waves there is no upper band edge, and for frequencies high enough, the wave 
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will be delocalized (ignoring Landau damping for the moment). 


It is not straightforward to compare parameters in the solar wind plasma with the critical ratio 
{W!V)c to determine dimensionality effects. The difficulty is a result of / in (4.2.3) being much 
smaller than the scale length of the change in i.e., much smaller than a typical wavelength. But 
on this scale the density fluctuations may not be (and in the solar wind plasma are not) decorre- 

lated, so that co^(jr = jl) is not independent from one lattice site to the next. Therefore, this situa- 
tion is significantly different from the solid state case where perturbations are (assumed to be) coin- 
cident with lattice sites. 


We can approximate the situation by imagining that our plasma is divided into a lattice with 
pieces of size where is the correlation length of the density fluctuations. Each piece is then a 
separate harmonic oscillator. 




f\/~^ 


\ 



X '2^4 

For a state with frequency o>, the wave number in each piece will be a function of the . The 


coupling between, say, oscillator 1 and 3 through 2 will be -p where is the group velocity for 

'c 


piece 2. Therefore, 


WIV = 




We use Vj rather than v when looking at the plasma on a scale I >X. When observed at 
small scales / the weak coupling v leads to propagation at the group velocity for distances 
>X. Further rescaling will not change the group velocity. Dividing the plasma into pieces larger 


than / - will smooth out the density fluctuations. On the other hand, on these larger scale 
lengths, the coupling between pieces becomes a complicated function to account for scattering and 
decreased transmission which are strongly dependent on frequency at this scale length. 
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Using an average group velocity given by 
( 





for Langmuir waves 
for electromagnetic waves 


and using the parameters for the 11 March event, estimating 1^ - 100 km. and — — = .01, for the 

Langmuir waves: 
k = 2.3x 10-* cm-*; 

Vj = 2.4 xlO’ cm/sec 
W!V=:\Qp 

for the electromagnetic waves: 

* = 1.7x10-’ cm-' 

Vj = 1.9 X 10® cm/sec 
W/V=1800 


In both cases, this greatly exceeds even the largest estimates of {WIV)^- Our estimates for WIV 

here may be too large because of overestimating and l} could conceivably be smaUer by as much 
as a factor of 25 or greater. 

Once establishing that WIV >{W!V)^ so that waves are strongly localized, of more practical 
importance will be estimates of WiE and £, the parameters used in the Escande and Souillard 
(1985) study. 

For the Langmuir waves 


WIE = 




= 1.4 


for the March 11 event where Vj, is the beam velocity. Likewise, if we could treat the transverse 
waves 1 -dimensionally (which we cannot) we would get, using Aj = 360 km for the March 11 


event: 
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fl] 

^ J Transverse 

If /j. = 100 km, E = 3.00, and plot 4.1 gives approximately y= .2. So | = 500 km. 

A difficulty in making localization estimates for the type III problem comes from the nature 
of the turbulent density spectrum, as discussed in section 2, which is that there are large amplitude 
long-wavelength fluctuations simultaneous with smaller amplitude, shorter wavelength fluctuations. 
If the localization length from the laner is smaller than the correlation length of the former, then it 
is the small amplitude short wavelength fluctuations that are important. That this may be so for the 
Langmuir waves is evident from the following. Using the 11 March 1979 parameters, if density 
fluctuations have a scale length 1^ = 2 km, then £ = 21. If the localization length ^~70 km, then 

we estimate from Figure 4.2 that WIE ~ .5, which gives — = .003, required to give this localiza- 

ft 

tion length. The actual localization is therefore dependent upon the fluctuation spectrum. 


Scaling Theory of Localization 

The first successful scaling theory of localization was proposed by Abrahams et aL (1979) 
who used the scaling of the conductance of a sample with scaling of size to find that all one- and 
two-dimensional waves are localized in a random system and that a mobility edge exists in these 
dimensions, delimiting a transition between extended and localized states. This scaling theory was 
verified numerically by McKinnon (1985), discussed earlier in this section. Although all states in 
one- and two-dimensions are localized, Figure 4.2 and the numerical work of Yoshino and Okazaki 
(1977) indicate that there is a transition between weakly and strongly localized states. 

Another way of viewing the localization problem mathematically is to look at the Green’s 
functions which are solutions of 


[-c^ -i- v(r) - £] G (r, Tq) = d{r - r ^) . 

They may be written as 


OA/R(r-ro,E)='^ 

a 


¥a(r) Va(ro) 
£ - ± 17 ] 
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where AIR determines whether it is the advanced (retarded) Green’s function, and a is an eigen- 
state label. Now G (r - ro) -» 0 as I r - tq I -> <» if all the wave functions are localized. 

As discussed in the previous section, the quantity 

D -<G'‘(r-ro.£)G*(r-ro.£)> 

where the average is an average over the statistical ensemble is important in the study of 

wave localization. For localized states 

-lr-rol/$ 

where d>\ is the dimensionality. An important property of this Green’s function is that as the 
disorder, WIE, approaches a critical value the correlation length | becomes infinite. This sudden 
transition of the correlation length is a universal property of phase transitions, in this case called 
the Anderson transition. By applying techniques IBrom renormalization group theory and the theory 
of critical phenomenon Wegner (1982, and references contained therein) argued that the rate of 
diffusion for frequencies above the mobility edge scales as 

£- [(£-£,)/£,j'' 

where v is some exponent (called the critical exponent). In 3-D one gets a “phase diagram” 
{.Stephen, 1983) 



where states lying within the hatched area are extended. In a solid, the boundary is determined by 
the dashed line, but as we have discussed, for a plasma, there is no upper band edge, and for 
sufficiently large values of ElV states will be extended. 


An important conclusion from this is that in 3-D, for any disorder there is some energy 
E^><o~-Sa>^, i.e., above the band edge, below which the states are localized. To see how this 
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might be the case we return to Anderson’s argument and consider the probability of convergence of 
the series 4.2.6. As i's-» £ + it] , where i} 0 and E — > -WI2, the probability of small denomina- 
tors decreases and therefore the probability of convergence increases. 

We consider the implications of this for the type in problem. The difference between the 
frequencies for the waves in the 11 March case are 

/t-/pi«».s45 Hz 
and 

/r-/pi»m .2 25 Hz. 

Therefore, even if we greatly overestimated W/V for the transverse waves for typical solar wind 
parameters, fr may be below the mobility edge in any case, since fj is so near the '*band edge/’ 
In the second place, if fr were in fact above the mobility edge, it would be expected that the rate 
of diffusion would be very small. The consequence of this is that a transverse wave that originated 
in a given region would still be “approximately” quantized, resembling a given stable bound state 
in quantum theory, and therefore, the implications for the type in problem would be little different 
than what is discussed in section 7. 

Localization and Mobile Density Fluctuation 

By applying normalization group techniques Vollhardt and Wolfe (1980, 1982) argued that 
the ladder and maximally crossed diagrams (4.1.7) are responsible for the Anderson transition. If 
these diagrams are the most important for localization then this provides an estimate for the effect 
of the mobility of density fluctuations. The destruction of enhanced backscanering in the weak 
localization limit was studied by Golubentsev (1983) who concluded that the loss of phase coher- 
ence due to the destruction of time reversal symmetry for the backscattered wave is important when 

T>T^ where t is the mean time between wave-scattering and is a measure of phase 

coherence time when the scatterers are mobile with the average time for a scatterer to travel a 
distance A, the wavelength of the wave. In this computation Golubentsev assumed that the scatter- 
ers had independent random velocities with a Maxwellian distribution. We estimate r = ^IVg and 
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Tx_ = and define 



where is the velocity of a density fluctuation. Then we expect localization will be preserved for 
» 1, and destroyed for ccq < I, 

If we use the ion-acoustic velocity for then for the 4 March event, for the Langmuir 
waves ao = .9 and for the transverse electromagnetic wave Oq - 100. This would indicate that the 
Langmuir waves may not be localized, but that the transverse waves probably are. We mention the 
following points: 

(1) As we shall see in numerical simulations in section 5.2., when localization is very strong, 
waves will become localized by a co-moving group of fluctuations and scattered by the remainder. 

(2) Although little is known about the source of the turbulence in the solar wind, very often 
in plasma turbulence there is an average drift along or across the magnetic field relative to which 
other turbulent motion is slow. Therefore phase coherence may be preserved longer since the ran- 
dom component of the velocity can be much less than the average velocity. 

(3) This estimate is based on only the destruction of time reversal symmetry of the enhanced 
backscattering. In fact, for strong localization, the phase reinforcement comes from a fortuitous 
coincidence in the scattering off of different fluctuations in the neighborhood of the localization 
peak. 

Density of States 

Also of importance to us will be the density of states. In this problem we are primarily 
interested in waves which are resonant with a particle beam. For this to be the case, we expect that 
the waves of interest will have a localization length that extends over many wavelengths. The den- 
sity of states will determines how many states/volume the beam can resonate with and therefore 
how much beam energy will be lost in exciting localized waves. 
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For any finite homogeneous system of length L, the density of states is given by the usual 
formula 


d^n 

dk^ 


2n 


1 d 


(4.2.7) 


where d is the dimension of the system. For a disordered system and for £ < £^, the density of 
states corresponding to bound states is discrete while for £>£^ there is a continuum of states 
corresponding to the freely propagating waves. The question is how closely we can approximate 
the density of states by the expression 4.2.7. As numerical work by Dean (1961) has shown, away 
firom the band edge, the density of states will closely resemble that of 4.2.7, while for frequencies 
near the band edge the density of states is highly perturbed. We expect that for waves whose local- 
ization length is many wavelengths in extent, the firequency will not be near the lower band edge 
(since the probability that the minimums of that density will correlate with an extended wave is 
very small). Therefore, we expect 4.2.7 will be close to the actual density of states. 

Using the parameters for the March 11 event, the number of states in one dimension in a 

Av^ 

localization length ^bIO km, resonant with the beam with = .05, the apparent width of the 

positive slope region, is approximately 




= 1.25 , 


2n Vj, 

Finally we consider the conditions under which localized states are resolvable. The density 
of states in 1-D is given by 


dn ^ 

dk~ 2n’ 


(4.2.8) 


2;r 


which gives dk = — dn. In order to determine the separation in frequency of modes contained 
within a localization length we start from: 


ct)“ = ci>^ 4* then 

day = dk where is the group velocity. Using 4.2.8 gives 
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Ao>= — ^ — . The length of time to resolve ^vo such states is given by the requirement 
Ao> >2;r, which gives which might have been guessed heuristically. 

Summary 

We have reviewed Anderson localization theory which leads us to expect that electromagnetic 
and possibly Langmuir waves wUl be localized by the density fluctuations in the solar wind plasma. 
Theoretical predictions are consistent with experimental observations. The theory in this section 
will be used in interpreting the numerical simulations of the next section. 



5. Numerical Simulations — Propagation of Gaussian 
Wave Packets in Random Media 


42 


We present the results of a series of numerical simulations of the time development in a ran- 
dom medium of an initially Gaussian wave packet. While numerical work has been done by others, 
in order to study localization the usual approach has been to diagonalize the Hamiltonian to hnd the 
eigenstates of a random lattice. This has usually been in the context of solid state theory and lat- 
tice sizes have usually been small (100x100 grid points) (Dean, 1961; Yoshino and Okazaki, 
1977). To our knowledge time development with mobile fluctuations has not been studied numeri- 
cally. 

Here we investigate the time development of an initial wave which presumably is the overlap 
of many localized eigenstates. In section 5.2 we also investigate the efifect of mobile density 
fluctuations. This study has used the program PERLOC, discussed in detail in Appendix B. This 
program propagates an initial Gaussian wave packet according to the wave equation 


dry/ X z 


2^ 

0 

with the initial conditions: 


^(:c./=0) = 0 

ot 

(/(jf.O) = v^ocos ko(x - jTo) 
periodic boundary conditions 

V/(0,O = v(Ar„.x,0. 
and where 

a>^(-v) is a random variable on the interval [ \ - dn , i ^ dn]. 

In most of the simulations shown dn = . 01, corresponding to a 1% density variation. The simula- 
tions are done on a spatial grid of 2048 or 4096 grid points. The fluctuations in are constant on 
steps typically of 4 or 8 grid points. A typical plot of the value of cof(x) is shown in Figure B.l. 
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We plot the following quantities at various limes, 

Psi: v^Cv, r) = amplitude. Since this is a linear problem, the units of this quantity are arbitrary. 


Energy: £(.r,r)s j + j 


dv(x,t) ] 

H 

d\y(x,t) 

a, J 

2 

dx 

L J 


This is the numerical value of the Hamiltonian density which generates the wave 
equation. This should be a slowly varying quantity, and therefore a better indica- 
tor of localization than y/ or individually, which oscillates in time as 

The following symbols are used in describing each simulation. The quantities are discussed in 
section 4.2 and values of W/E and E are used to enter Figure 4.1 to determine the predicted locali- 
zation length 


central wave number of the initial Gaussian wave packet 

c = velocity in the dispersion relation a>“ = + c~k“ ; co“ s 1 

He = number of grid points in a density fluctuation step 
= length of density fluctuation step (correlation length) 

£ s k§ 1} 

WfE s Scopk^c^ 

y= Liapunov exponent from plot 4.1 

^ = /^/7= predicted localization length from the theory of Escande and Soiullard (1984) 

dn: the plasma frequency is a random variable on the interval [I - dn, I + dn] 

kc^ 

= (group velocity) 

Vgt= total distance propagated at the group velocity at the end of the run. 

In interpreting these simulations the following general remarks pertain. 

1) If then energy - Therefore, the localization length will be the '/z width 

of the wave packet at the amplitude E = £^„ 

2) Although it might be expected that results would only be valid after a time t such that 

in fact we find very good agreement with theory when v^t 

3) An important qualitative difference between a freely propagating and a localized wave can 
be seen in the energy. For a propagating wave, y/ cos(kx - co^t), and so, for small 

E — ^ [cos‘(^-co^O + sin‘(A*.v-a)^/)] - constant. For a standing wave, on the other 
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hand, iff - cos kx cos cot and so for small E co^ cos“ kx(cos~ cot 4 * sin“ cot)- cos“ kx. For 
a propagating wave the energy density will have a smooth envelope, while for a standing 
wave it will be modulated with wave number k. 

Free Propagation 

Figures 5.1 and 5.2 show freely propagating waves (constant density). These are to be com- 
pared with the other simulations in this section which are all with randomly fluctuating densities. 

These two runs show the free propagation at the group velocity of two initial Gaussian wave 
packets. The initial excitation of the form 

¥(x) = ¥o cos ko (x - XqY, \j/(x) = 0 

is seen to split in two and propagate in opposite directions. The dispersion relation is 

k^ with (of=i 

in both cases Vgt is very close to the distance which the p)eak of the envelope propagates. The 
parameters in these simulations are the following. 

Figures 5.1.1— 5>1.5 Figures 5>2.1— 5.2.5 


A-o=.16 

*0=1-7 

II 

c=.05 

V =.0064 

v^ = .0043 

v,/ = 51.2 

Vjj/=34.4 

/I. =4096 

n«=4096 

II 

dt = .2 


5.1 • Localization Length Studies 

The following series of simulations shows the effect of time-stationary random density fluctua- 
tions for a wide set of parameters. The observed localization lengths are consistent with the values 


from Figure 4.1. 



I 
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Figures 5.3.1— S.3.6 
Parameters: 


This set shows the time development for the same Gaussian 
wave packet as for the freely propagating case shown in plots 5.2.1— 
5.2.5. In this case, however, the density fluctuations are strong and 
localization is quite evident. In plot 5.3.2 we see the high wave 
number small amplitude signal propagating at the velocity c. This is 
the Sommerfeld precursor {Jackson, p. 9, p.313 £f.). With WIE =9.8 
this is far removed from the plotted region. We might estimate y- .1 
and therefore - 12, which is of the same order as the observed localization length. In this case, 
however, we do not expect k^lc to be a meaningful quantity since the wave is so strongly distorted. 


)to=.16 

c = .2 
dn = .01 

/If = 16 

/f=1.22 

£ = .04 
W/£=9.8 
v, = .Q064 
/I, =4096 
d, = .2 


Figures 5.4,1— S.4.4 

Parameters: 

ko=.32 

f = .2 
d/i = .01 

;if= 16 

/f = 2.54 
£ = .66 
WIE=2M 
v^ = .0128 
v'r=307 
y=8xl0~- 
^=31 


Figures 5.5.1— 5.5.6 

Parameters: 

Ao=.85 
c = .l5 
dn = .0l 
/If = 3 
/f = .46 
£ = .154 
W/£ = .61 
y=3xl0“^ 

^ = 153 
v^ = .019 
v,/=570 


We have run this for 120,000 time steps to f = 24,000. Despite 
the very strong localization, there is still very good agreement 
between the predicted and the observed localization lengths, y is 
estimated since W/E = 2.44 is not on the graph of Liapunov 
exponents Note that localization is very strong even though 
v ,//4 = 10 . 

» 


We show the time development of this in some detail. This 
shows a case where the final localized wave packet is -40 
wavelengths long. The initial wave packet is confined to a region 
much smaller than the final localized wave. The wave packet splits 
and begins to propagate, as seen at / =2,000, however, there is evi- 
dently already significant reflection and a significant amount of wave 
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energy remains in the initial region. By t = 12,000, the ounvard pro- 
pagating portion has been significantly depleted and by time 
t = 18,000 there seems to be little propagation. It is evident in the energy plots at t =20,000 and 
r = 22,000 that the observed localization length is -150, consistent with predictions at t =28,000. 
It is more confused since propagated energy does wrap around due to the periodic boundary condi- 
tions. 

The initial and “final” conditions are exhibited for the following sets of parameters. 


Figures 5.6.1— 5.6.3 Figures 5.7.1— 5.7.3 


Parameters: 

Parameters: 

A - o = 1.28 

A-o= 1.7 

c = .07 

c = .05 

dn- .01 

dn = . 01 

rt^ = 8 

n^ = 8 

l,=.306 

/^=.306 

£ = .15 

E = .21 

W/E = l.2A 

WIE = 137 

v^ = .0062 

v=.m 

Vgt = 50 

v,/ = 34 

r=8xl0“^ 

r=i.5x 10-- 

1 = 38 

1=20 

Figures 5.8.1— 5.8.3 

Figures 5.9.1— 5.9.3 

Parameters: 

Parameters: 


-to=l-7 


c=.05 

dn = . 01 

dn = . 01 

w^=8 

«<-=8 

/,= 1.23 

4 = .3 

v, = .013 

v, = .004 

vy=230 

v,/=34 

£=2.47 

£ = .27 

VF /£=.61 

WIE = \3i 

y=1.2x 10“^ 

r=2x 10"- 

4 = 100 

4 = 15.3 
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5,2. Effect of Time Varying Density Fluctuations on Localization 

In a real plasma the density turbulence is not time stationary, and therefore the effect of 
motion of the density fluctuations on localization must be studied. As discussed in section 4.3 we 
expect that the parameter 




^2 c2 


1/3 


will be important in studying the effect of the mobility of density fluctuations. 


One would expect that for :s> 1 waves will remain localized, but might spread at a rate 
corresponding to the velocity spread of the density fluctuations. As decreases to ~ 1 one 
might expect to reach a point at which the wave is no longer localized and therefore the wave 
packet will spread diffusively. In fact, in the numerical work presented, for strongly localized 
waves, this does not seem to be the case. 


Density fluctuations moving uniformly with a velocity v< v^, the group velocity of the waves, 
are stationary in a co-moving frame and should not affect localization. This has been verified 
numerically, but is not presented here. 

A second possibility for one-dimensional turbulence is that density fluctuations have a con- 
stant speed and may be directed either to the left or to the right. To investigate the effect of this 
we form two random independent distributions 8 n^{x) and on the interval [.5-dn/2, 

.5 + dn/2] and form the time dependent distribution 


8 n (xj) = 8 n^ {x- c\t) + 8 /I 2 U + cj ) . 

The results of these simulations are shown in Figures 5.10.1—5.10.6 and 5.11.1—5.11.6. In these 
cases one has to remember that the sum of two uniform distributions is not a uniform distribution, 
and so the results of previous localization length predictions do not apply directly. 

For these we use parameters as before except that: 


— ^ = -( 0 - + 3c, 


dr 


dx- 
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so V, = 


3k c: 


Figures S.10.1— 5.10.6 
Parameters: 
<//t]=<//>2 = .004 


y-.l 
^ = 10 
ac=1.2 

c=.04 

c,=.004 

Figures S.11.1— 

Parameters: 

ito=.426 
c=.2 
dn = .05 


In this case we see that the initial wave packet splits and 


II II 

separates at the 

v.=.003 

v,/=42 

moving density 

vq p 
II II 

scattering by the 


In this set we show the resulting distribution of £ in the cases of 
oppositely moving density profiles with: 


«,=35 
/^= 10.7 
v,=.05 
v,/ = 1635 
£ = 18.5 
W'/£ = 1.6 
y= 

1 = 50 

(from Fig. 5.11.2) 


Figure 5.11.1—5.11.2 

Cj = 0 (stationary profiles) 

ac = oo 

Figure 5.11.3—5.11.4 
c, = .005c {cjv^ = 1/50) 
(Xq = 8.7 

Figure 5.11.5—5.11.6 
c, = .01c (c/v^ = 1/25) 
oc = 5.4 


In these cases the wave packet spreads at the speed c^. Most of the wave energy seems to 
remain confined within a space determined by Xq±(^ + c^t), although we do see indication of some 
additional spread of wave energy, especially in Figure 5.11.6 in the regions jr< 185 and v>415. 


In the following simulations we assume density fluctuations with widths 1^, amplitudes, random 
on the interval [-dn,dn], and velocities random on the interval [-c, cj. In other words, each 
density fluctuation has an independent random velocity. 



Figures S.lla.l— 5,12c.l 
Parameters: 
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io=.426 
c = .2 
dn = . 02 
n^=14 
/^=4.29 
v^=.05 
v,r=480 
£ = 3.3 
W IE = .92 
7= .06 
4=75 


Figure 5.12a.l— 5.12a.2 

c, = 0 (stationary density fluctuations) 
Figure 5.12b,l-5.12b.2 
c, = .01c {cjv^ = 1/25) 
ao = 4.14 

Figure 5.12c.l-5.12c.2 
c, = .01c {cjv^ = 1/2.5) 

«G = *9 


In this case again the spreading of the wave packet seems limited by c^t . 


Figures S.13a,l— 5.13c.3 


Parameters: 

c = . 12 
dn = . 02^ 

/, = 2.45 
v^ = .037 
v^/ = ll80 
£=4.7 

7= .03 
<^ = 85 


Figure 5.13a.l— 5.13a.6 

In this sequence we show the time development for stationary 

density fluctuations. 

Figure 5.13b.l-5.13b.3 
As above, but c^ = .0025c (c,/v^ = 1/3333) 

= 6.8 

Figure 5.13c.l— 5.13c.3 
As above, but c^= .01c (cjv^ = 1/31) 

aQ = 2J Note that the vertical scale is different in 5.13c.l and that 
the localization region is the same as the others. 


Figure 5.13d.l— 5.13d.3. As above, but c,= .05c (c^/v^ = 1/16) = .93 


Figure 5.13e.l— 5.13e.3. As above, but c, = .Ic {cJVg = 1/3) aQ = .58 


Again, in this sequence the spreading of the wave packet seems to be limited by c^/, despite 
small values of order to decrease further, it would be necessary to increase Cj~ Vg, At 

this point localization ceases to have any meaning. 


Conclusions 

We have found excellent agreement between the localization length predictions of Escande 
and Souillard and the confinement of waves observed in numerical simulations. The interpretation 


of 
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ko s ^{(o^-(oj)fc^ 

as the central wave number of the confined wave packet is borne out, and does lead to, the correct 
localization length predictions. 

Strong localization seems to depend weakly on global properties of the medium and strongly 
upon local properties. By this it is meant that the structures and spectrum of a localized packet is 
largely determined by the properties of the medium where most of the wave energy is confined, and 
the effect of the medium far from this seems only to insure that what little wave energy does reach 
distant points is scattered back. It is this property which seems responsible for the fact that a local- 
ized packet, in the presence of mobile density fluctuations, spreads at the speed characteristic of 
the speed of the density fluctuations rather than at, for example, 

We see for example, in the early simulations and in 5.13a.2, that initial spreading occurs at 
the group velocity, despite the existence of density fluctuations up to the point where the wave is 
confined. 

It must be kept in mind that the results shown here for mobile density fluctuations are not 
necessarily valid for higher dimensions, and in particular it might be expected that waves near the 
mobility edge in three-dimensions are delocalized by time dependent density fluctuations. On the 
other hand, the apparent lack of dependence of localization in one dimension on the motion of den- 
sity fluctuations provides hope that things will not be radically different, and the argument is 
further supported by the weak dependence of localization length on dimensionality for strong local- 
ization as discussed in section 4. 
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6. Beam-Plasma Processes in Random Media 

In section 6.1 we develop a theory for the weak-beam plasma instability in random media. In 
section 6.2 the results of numerical simulations of this instability are presented. In section 6.3 we 
develop a theory for the saturation of the instability. 

6.1. Theory of Beam-Plasma Processes in a Turbulent Plasma 

If waves are localized in a turbulent plasma as theory and our simulations suggest, then a 
beam-plasma instability could be expected to excite localized waves with a central wave number 

CO 

resonant with the beam particles. 

In order to understand how this interaction comes about, we return to the Lagrangian discus- 
sion of section 3. We shall concern ourselves with the ponion of the Lagrangian 

IBO 

We make the assumption that we can represent the fields in this expression as a sum over 
localized waves whose representation is an exponentially localized wave packet, i.e., 

m 

where 

t) = (t>„ cos co„t cos k„(x-x„)e~’‘"' 
where k„ = 1/|„. 

As seen from the simulation in the previous section, this is only an approximate representa- 
tion; however, it will allow us to make some theoretical predictions, and it should be reasonable for 
waves where the localization length is much greater than the central wave length. 

We can rewrite this as 

f r 

J COS [ ^'(-V " ] + COS [ ^(.r - ,V„) + 0>„t] 


Pi Xi--Z <P(X,) 

Zm 
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where 


(k - k„f + K- 

From this we readily obtain the generator for the Lie transform 


W = 6i jdt 

i 

^-'Lei^jdkSJk) 

i.m 

^ sin [^( at, -. 0 + <o„r] 

kVi + to„ 


sin [k{Xi-x„)-a>„t\ 
kv: - 0)„ 


If ^ €i <t>(x,) where the sum is over all beams and plasma particles, then 


a-< 2)= I /'w, X ^:^fdk'dk"s„(k')SAk") 

2 2 J „ „ (27t)- ■' 

I sin [k' (.V, - x„) - co„] ^ sin + <>>„/] 

\ k'Vi-(o„ k'v,+co„ 

cos (* " (.V, - x„) - <o„t] + cos [A: ” {x, -x„) + <aj\ 

/ 

Evaluating the Lie bracket and taking a time average over the rapid phases to„/, coj 




K(-^=Y,-^-^\dk' dk'' S„(k’)S„{k'')k’ k'' 

tH, J 


sin k Xff^ sin k (Xj Xu^') sin k (,V| 'Vn,} sin k (x^ x^') 

(k'v, + co„f , (Ar'v,-‘‘>m)' 

Next we assume that, in this representation, we may replace the sum over particles by a spa- 
tially uniform average. Using the Vlasov representation: 


X = X d^vS^ix - Xi) 5^(v -\>i)= \d^x d^Y noS(y) ■ 

t i 

The assumption that we make here is that once we have taken into account the effect of den- 
sity fluctuations on the wave propagation which is to localize the waves in the problem, we can use 
average medium properties (such as averages of the dielectric function) to account for wave-wave 
and medium-wave coupling effects. The simulation results that will be presented in section 6.2 
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bear out this assumption. 
This gives 



2n m 




X \d^k' dH” k' k” S„ik')S„{k”) 


S{k' + k")-S{k' -k") 
2 


(k'v + o>„)^ (k'v-ioS- 




X Uv^(v) 5 

rn J ik'v + o>S- (k'v-o>„)^ 


Likewise we can write 


^fields = -^'L^mjdkk - Si (k) 


^ 8;r- 

We can then write our Lagrangian as 


L<2) = ie<^> + Z. 


fields 




npe‘- 

2n 


1^3, ^(V) 


g(v) 


(k'v + o)- (A-'v-co,)- 


J_ 

An 


=-S 


72 I 2 


Sir 


^ (f (^/nt ^ 

where e(^,co) is the usual plasma dielectric function for a plasma with uniform density rip and 
velocity distribution ^(v). 


But then the action is 


S<’-> = pt a,]. 

Varying the action with respect to <p„, and setting it to zero: 


0 = 


SS^^> 

5 


gives 
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[e (k„, co) + e(-k„, o>)] = 0. 

We note that with the beam present, 

eik„,o))^e(-k„, (o). 

Following the standard analysis this will give a growing wave mode provided 


dv 


^ 9v 




> 0 . 


For all the cases we are interested in the beam velocity 10 v,^. The damping from the 
Maxwellian is negligible. Therefore, (j>„ should grow at a rate the standard Landau growth rate. 

We therefore expect that in a turbulent plasma those localized modes which are resonant with 
the positive slope portion of the beam will grow with the same growth rate as a plane wave with 
the same central wave number in the homogeneous case. A major difference between the two will 
be that the localized waves will remove much less energy from the beam for any given maximum 
wave amplitude. 

For the sake of completeness, we can also look at single particle motion. If we have the new 


single particle Hamiltonian K^~\ then the guiding center motion of particles is F, = - 


^2 j ^ j 2 

-T S — dk"s„(k')s„(k")k'k" 

2 “ OT, ( 271 )- •' 


dX, ■ 


sin k (x, ^tn k 


1 


1 




This can be evaluated noting that: 


jdk”k" S„(k”) sin k"y = [cos k„y 

and near resonance for v,, 


jdk'k'SJk') 


1 

(^w ~ ^nt ''V)” 


sin k'y 
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7t 9 


9 


f \ 
0> 


V 

V ' / 


0> 






COS 


eg 

V; 


[©()-) -0(-y)] 


where 


] 1 if y>0 
[0 if y<0 


Far from resonance when 


a>n - I k^Vi 1 






' . ' 


iKvi-(o„f (A-.V, + <•),)- 

This is just the usual ponderomotive potential (Schmidt, 1979, p. 51). 


Near a resonance. 


<^n - I ^‘,v, I < K, 


‘ 2 n,,(2n)^ 


7t 9 


->r^ cos k„ (k - jr,) 


9 


o> 


vS 




y 

V ) 


cos — (-V-.V,) 


|^©(.r-.v,)-0(.v,-.v)j 


8 m; vf 


K 

9 

/ 

( \ 


( ‘ 
COf, 


V; 


9 " 


( / 



k 



sin 




*/<«•+ 1 

j^0(.r-.r,)-©(jr,-.r)j| 




fr 



2 > '> 

ei <l>« 




Sm.vfx:, 

- k„)-/K; 

^ COS 

[ 






X [©Cv-.vj-0(.t,-.v)] 

We note the following points. Everything outside of the braces is dimensionless. AT/-' - — 

and therefore is strongly dependent on the localization length. It also changes sign with a>„ - k„ »•„ 
i.e., as the particle is faster, or slower than the resonant velocity. We will return to this subject 
when we discuss saturation mechanisms for the beam-turbulent plasma instability in section 6.3. 



56 


6.2. Simulation of Beam*Plasma Processes 

In this section we show the results of a series of simulations of beam-plasma instabilities. 
The method of simulation and the program are described fully in Appendix A, but briefly is as fol- 
lows: 


The beam is simulated by particles which obey the equations of motion: 



= — £-U,) 

m 


«be.m (x.r) = ^ 5 (x - .t,(r)) 


For the background, define 


with equation of motion 

-:^ = co;(.Y)i/ + 3vr^ 
dr ax“ 


where a>^(.t)s4;r;i, {Xyt)e“lmg is a random variable. 

For the fields: V ♦£ =47rp, where p = -e/Jbc^(jr,0+ 

We plot E and as before, except that now represents the background charge density per- 
turbations due to plasma oscillations. The plot parameters are as in section 5 with the following 
additions and exceptions. 

V,, = effective electron thermal velocity of background plasma 
njnp=^ beam density/effective background plasma density 
Vi, = beam drift velocity 
= beam thermal velocity 
A'o= beam resonant wave number 


V. = 


3A'o V,; 


= group velocity of beam resonant waves in the background plasma 


# of states = 


k^L 

2k 


( .. \ 


V. J 


, which represents the expected number of localized states which 


should be excited by the beam. (See the discussion in section 4.2 on density of states.) 
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The initial conditions for all of these simulations are va ( a ‘) = 0, \^(. v ) = 0, and beam particles 
uniformly distributed in space with average velocity and a Maxwellian distribution with thermal 
velocity \s. 


This set shows a system without density fluctuations and should 
be compared with 6.3.1— 6.3.6 which has 1% density fluctuations. 
The noteworthy points in this are the smooth variations in the 
envelope of y/, and the purity of the harmonic content. The envelope 
modulation can be understood as the beating of a narrow spectrum of 
waves excited by a beam of finite temperature. If we take = 2;r, where ^ is the characteristic 

length of the modulation and Ak = A ^o>/vj = co^ = .066. Therefore ^ = 94, consistent with the 
characteristic length of the modulation of yr- The other noteworthy point is that the energy density 
plot is characteristic of a traveling wave. Furthermore, comparing the positions of the peaks at 
t =3,000 and t =4,800, it is apparent that they have moved a distance VgAt = 18 between these two 
limes, as one would predict. 

Figures 6>2,1— 6.2.15 

6.2.1 shows the distribution of density fluctuations. 

6.2.2 shows the initial beam-particle phase space (velocity vs. posi- 
tion). 

6.2.3 shows the initial beam distribution function. 

At r=0, v/ = 0, and Figures 6.2.4—6.2.15 show the subsequent 
time development of the beam-plasma instability. 

At t = 200, y/ is approximately uniform in space, and the energy 
shows that the wave is a traveling wave. By t = 1200 the localized 
structure is beginning to show. Between times i - 4,000 and 


Parameters: 

v,=.035 

Vft = 1.0 

rii,/npi = 9xl0-^ 
dn = . 01 
nj.=4 
)to=1.0 

/, = .6l 
£ = .38 
W IE = 2.11 
y=.06 
^ = 8.7 

v, = .0036 
v,r = 20 
# states =5 


Figures 6.1.1— 6.1.3 

Parameters: 

v,=.07 
Vft = 1.5 
v„=.15 

V"^/=9xl0“* 

ko=.61 
v. = .01 
dn=0 
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t = 5,000, the mode structure is reasonably stable and the growth has 
saturated. 


Simulation shows the same initial conditions as in 6.1.1— 6.1.3 
except that 1% density fluctuations were introduced. The energy den- 
sity plots are consistent with standing waves. 


Figures 6,3.1^>3.6 

Parameters: 
vv=.07 
V£, = 1.5 
v,,= .15 

ni,/nf,, =9x10-^ 
4/n=.01 
ko=.61 
ric = S 

/<,=2.34 
£ = 2.46 
W/£ = 1.5 
y=6xl0'^ 
^=27 
v, = .01 
v^t=58 
# states =13 

Figure 6.4.1— 6.4.6 
Parameters: 

vj, = 1.26 
v„=.12 

ni,/npi = 9xl0~^ 
dn = .02 
ko=.19 

«j.=4 

/^=i.l7 
£ = .854 
W/E = 1.06 
r=2xl0"^ 

4 = 59 

v^ = .024 
v,/ = 142 
# of modes= 15 


Figure 6.5.1-6>5,2 
Parameters: 
v,,=.04 
Vi =.7 
^=•1 
dn = , 01 

nj rip, =9xi0~^ 

1:0=1.43 
"c=4 
/,=.6 
E = .ll 
WIE = \.Q 
7=2x10-2 
4=30 
# states=10 
v, = .006 
V=26 


Figure 6.6.1— 6.6.2 
Parameters: 
v,=.05 
v* = .7 
v,=.08 
</n = .01 
/ii/«P=9xlO-* 
ko= 1.43 
«c=4 
/^=.61 
E = .ll 
WIE = 1.9 
7=5x10-2 
4 = 12.2 
# states =8 

v, = .0036 
v,/ = 13.6 


Figure 6.7.1— 6.7.2 
Parameters: 
v,,=.15 
Vi = 1.5 
^=•1 
</n = .01 
''fr/'»p=9xl0-* 
ko=.66 
n,=4 
l,=.6l 
£ = .16 
W/E=2.3 
7=6x10-2 
4 = 10 
# states =2 

v^ = .0067 
v,/ = 13.2 


Growth rates and levels at saturation 

The peak energy density (£) versus time is plotted in Figure 6.1 1 for the simulation presented 
in Figures 6.2.1-6.2.15, giving an observed growth rate between /=200 and / = 1200 as 
y/o>^=.00l5. The growth rate between / = 1200 and /=3500 is y/ = .00024. (Note that the 
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growth rate for E is 2y.) The standard expressioa for the Landau growth rate of a weak beam is 


y _ 1 ^fb 






tolk 


{Cheriy 1976), and gives for the fastest growing wave resonant with a Maxwellian beam: 


Tmax 

<o„ 


S.38 


A. 

co; vm. 


= .00034 


for the fastest growing wave for the parameters of this simulation. We had predicted that the 
observed growth rate for a localized wave should be y/2=,000l7. We note that the observed 
growth rale between r = 1200 and / = 3500 lies between y/2 and y. That the observed growth is 
larger than y/2 is to be expected since the growth time l/y--^/v^ and so the resonant wave is only 
partially reflected during the growth time. The very large growth rate for ^<1200 is a result of 
discrete particle effects early on and the “quiet start” beam initialization {Birdsall and Langdon, 
1985). Resolving these difficulties would require many more particles and a much slower growth 
rate, straining available computer time. 

We can scale the levels of saturation in our simulation to values for the type III solar radio 
burst in mits units as follows. Let primed units refer to type III parameters in mks units and 

unprimed quantities to simulation values. Then from the equation of motion x-(qlm)E the scal- 
ing is 


E' im'/q') j/ ± 

E (rnlq) v t' 

We use ///' - cop/o)^, rn'lq' = electron charge/mass ratio in mks units and 

v' _ Type III beam velocity (a/?/:.?) 

V Simulation beam velocity 

Since E -Anpik = x^fik; k = 1; and at saturation (Figure 6.2.12) y/-2x 10"*^, the saturation field is 
£' = 3 mV/m 

which should be compared with the 1 mV/m observed field for the March 1 1 event. This scaling, it 
should be noted, is only valid for the case that and V/,/vv>,^ are the same for both 
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plasmas. This is not exactly the case, and furthermore, as discussed in the summary to this section, 
the distribution function for the solar wind is not a Maxwellian. 

6.3. Saturation of the Beam-Plasma Instability 

A fundamental problem in understanding the type III solar burst is the mechanism for satura- 
tion of the beam-plasma instability, limiting growth to the observed levels, despite the persistence 
of the observed positive slope at a distance of 1 AU. An extensive study of this problem was done 
using quasi-linear theory (Magelssen^ 1976; Magelssen and Smith, 1977). These authors concluded 
that an important effect in the persistence of the beam is that at a given point in space, the beam 
velocity decreases with time. Slower beam electrons later in time reabsorb plasma waves which 
were resonant with the faster electrons that produced them earlier in time. At 1 AU the ISEE-3 
data do not support the quasi-linear relaxation hypothesis, since a strong positive slope is observed 
on the electron distribution function. 

Papadopolous (1974) has proposed that the oscillating two-stream instability is responsible 
for removing resonant Langmuir waves, thereby limiting the growth rate of the beam-plasma insta- 
bility. Observed wave levels at 1 AU are below the threshold for this process, however {Lin et aLy 
1986). As discussed in section 4.1, Muschietti (1985) examined the effect of Langmuir waves 
being scattered out of resonance with the beam electrons by solar wind density fluctuations and 
concluded that the process was “too efficient.” They found that unless special assumptions are 
made about the anisotropy of the density fluctuations the beam plasma instability will be 
suppressed. 

In the Whelan and Stenzel experiment (1984), the beam was observed to travel the length of 
the experiment with no evidence of quasi-linear relaxation, although Langmuir waves from the 
beam-plasma instability were observed only in the first few centimeters of the experiment. They 
concluded that particle trapping was responsible for saturating the instability. In their simulation of 
this problem, Pritchett and Dawson also concluded that trapping was responsible for saturation of 
the instability. In both cases, however, the beam was relatively cold. 
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Davidson (1972, p. 55 ff.) gives the characteristic trapping time for a particle resonant with a 
Langmuir wave with electric field E and wave vector k as 

for the 11 March event this gives 

r^ = .(X)16 sec . (6.3.1) 

On the other hand, the transit time through a localized wave packet of 70 km extent using from 
the 11 March event is t -^/Vh = .002 sec. It should be kept in mind that in 6.3.1 we have used the 
peak field. Not only is the average trapping time much larger due to the lower average field, but 
the bounce frequency changes radically on transit through the localized structure. Therefore, parti- 
cle trapping does not occur, and we propose another saturation mechanism in what follows. 

We show the distribution functions at saturation for two cases. The runs are similar to those 
of runs 6.1 through 6.7, however, we have used a short simulation volume (L = 150 vs. L = 1200) 
in order to increase the phase space density of simulation particles. This results in much smoother 
distribution function. 


Figures 6.8.1^.8.3 
Parameters; 

vv=.15 

v,=2.2 

v„=.15 

V«p=9xl0-® 
dn = .02 
n, = l2 
y=.06 
^=30 
ko=A5 ■ 

/^ = 1.8 
£ = .65 
W/E = 1.5 


Figures 6.9.1 — 6.9.3 
Parameters; 

v,=.13 

Vfc=1.8 

v,=.15 
= 9x10-* 
dn — .03 
«,= 12 
y=.09 
^ = 20 
l:o=.55 
/^=1.8 
£ = 1.0 
W/£=1.8 


The important result is that in both cases the distribution function is fiattened in the vicinity of 
the wave-resonant velocity. There appears to be no overall healing or diffusion of the beam, a 
characteristic shared with every simulation that we have done. 
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The crucial difference between our condition and quasi-Iinear diffusion is that here the parti- 
cles see only one wave at a lime. Furthermore, the localized wave, although having a width in k 
space is not well approximated by the random phase approximation, since, in fact, it is a coherent 
structure, the phases are well determined, and first order changes to the particle velocity are enough 
to maintain a certain phase relationship between the wave and a particle. 


Velocity space diffusion of particles by localized wave packets 

In order to investigate the effect of a localized wave packet on the distribution function, a sim- 
ple program was run. We followed the orbits of an initial particle distribution governed by the 
equation of motion 


— - = Eq cos / cos .r e~ 
dt“ 


where k: is a variable parameter. was chosen so that the trapping velocity 



meaning that a particle of velocity Vq± .1 would be trapped by the wave if k = 0. The results were 
as follows: 


Figure 6.10.2 the initial distribution function for the numerical experiment was: 

( constant .8 < v < 1 
0 otherwise 

Figure 6.10.2 distribution function after passing through the wave packet for v = .1. 

Figure 6.10.3 as above but fc = .2 
Figure 6.10.4 as above but k = A 

It is clear that the particles have diffused in velocity space around the resonant velocity 
Vres= 10. The diffusion is significant in the neighborhood of velocity v= 1.0 with width approxi- 
mately = ± k:. In the case (not shown) where 
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constant .8 < v < 1.2 
0 otherwise ’ 


no distortion of the distribution after passage through the wave was noted. 
In order to explain this we begin with the Vlasov equation: 


We expand / and H by order in the field: 


(6.3.2) 


/=/(0) + /0) + /(2) 




Then 6.3.2 becomes 


^ = or, 

/^’^ = and 






where the bar indicates time-averaging. 


Using our previous expression for w, 


= -^\dk [ ^ 


sin {kXi + or) sin (/br, + or ) 
kVj + o k\’i - o 




jdk' dk" 


[(ko-k')^K^-] [(k^-k"y-+K'-] 


„ I C0S(A-".V + £ 0 /) , cos{k"x-ot) 5/o 

X 1 rrr~ + — tt, — . cos k .x cos or 

I k v + o k V- o dp 


_ -e-(p- (ji^> K-k'k" sin ^-'.v cos ;-".v ^ 

A ir~ - f.. ...^ . ... .1 


{ko-k')- + K^- r(A-o-A-'')- + K-l 


1 _ _1 ^ o_ 

k"v^(o k"v - CO \ f)p 
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Keeping resonant terms and ignoring ponderomotive effects, this can be evaluated as 

• J ^ -KTltl 

/ X fcv Sin Ato^v Sin — x e 

ekp4> ' d_ I *' I 3/o 

dt m J 9v (o) - *ov)^ + KT^v^ 3v 
This is the standard form for a diffusion equation. Although we might try to integrate this 
numerically, it is instructive to examine this expression and make some estimates of the diffusion 
rate. We note first that D(v) will be significantly different from zero only in a region 

<v < 
and 






As an example we use the 11 March 1979 parameters. For v = .lc, k- = 1/7x10‘*iii (20 
wavelengihsAocalization length), and /co = 2x 10“^/cm which are values typical of a xypc III burst, 
we obtain = .5 mV/m. This number should not be taken too seriously, but is close to the 
observed - 1 mV/m of the 11 March 79 event. 

Likewise for our simulations (Figures 6.2.1-6.2.15), equation 6.3.3 with ^*£ = v^ gives us 

= v^. Using the simulation parameters m/e= 1, jc = 1/<| = 1/8.7, v= .1, we predict at satura- 

e 

tion i^= 1.3 X 10“^ which is within a factor of 2 of the values in Figure 6.2.12. 

This theory for resonant diffusion is similar to quasi-linear theory {Davidson, 1972, p. 151, 
ff.), except that we do not make the assumption is usual in weak tur- 

bulence theory, but keep the Fourier representation of the waves in the system and integrate over 
that. 

Summary 

The consequences of wave localization in one dimension as presented in this section are con- 
sistent with satellite observations of the weak-beam plasma instability in the solar wind. Two con- 
siderations have been ignored in this section: 

i) The solar wind is three-dimensional. This is discussed further in the conclusions (section 

8 ). 

ii) The solar wind plasma has a non-ihermal tail (see figure 2.3). Among other things this 
reduces the positive slope of the distribution function reducing the growth rate. The 
growth time 1/y, in our simulation are < the localization times, while growth times in the 
solar wind plasma are ~l sec. 

For the 11 March event v^ = 2.4x 10^ cm/sec. If -70 km, the localization time -r = .3 sec. 
So yr = .3. If anything, this should mean that localization is even more evident in the solar wind 


plasma. 
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Simulating the non-ihermal tail on the background distribution function is a ver>' dilficult prob- 
lem because of limitations on the number of simulation particles (see Appendix A). The experi- 
ment done by Whelan and Stenzel (1985) suggests that this non-thermai tail may result from a 
strong turbulence mechanism, which would have significant implications for the mechanisms dis- 
cussed here. These mechanisms may be important far upstream in or near the solar corona, but as 
discussed in section 2, are inconsistent with satellite observations at 1 AU. 
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?• Three-Wave Processes 

Although this section is last in the thesis, the considerations in it were the motivation for 
studying everything so far presented. In section 7.1 we will review the theor>' of 3-wave processes 
in homogeneous media, and then in section 7.2 discuss the convective problem as it applies to the 
type ni problem. We conclude that wave convection in homogeneous media inhibits the 
I-^r + 5 decay process. We review previous work on 3-wave processes in inhomogeneous 
media. In section 7.3. we develop theory for 3-wave processes when wave localization occurs. 

7.1. Three-Wave Decay Processes 

In the absence of density fluctuations, the second 
written as 

where the first integral on the right is of second order in 
Starting from the expansion (3.4) 

= 1 z. ^ //(O) + 1 £ 2 ^(i) + £,//(2) 

and using 
we obtain: 

If we can use the eikonal approximation to e.xpress the potentials 
0=^0^. + cc (complex conjugate) 

k 

and 


and third order terms in the action can be 



the fields and the second is of third order. 
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then, as Kaufman has shown 

5^-^ = -^ ^d^xdt e^{k,(o\xd)^\ U,0+ X ^ (k,co,xt) e^,(k,x) 


Ab(xj) 


where £ is the dielectric function for electrostatic waves, D is the electromagnetic dispersion ten- 
sor, and if, is the polarization vector of the electromagnetic wave, b. 

We have made this separation because we shall be concerned with electrostatic and transverse 
electromagnetic waves. As Kaufman has shown, one can use this to find the variable canonically 
conjugate to the phase of the wave in the eikonal approximation localizing ^ttd 

m = — — , one finds 


4Cv,r) 


5S 


S[dBjdt] 


Likewise 




ss 


5S 

5(0 


8S 


de^ (xj;k^,(oJ 
d(o 


<t>a 


/ 

’ 11 

o 

3o)fc 

el D (x,r,k,co)et 

J > 


\Af, U,/)l' 


5[d9f,/dt] 5(Of, 

These expressions will be recognized as the standard action densities with the Hamiltonian 

Ct 

We may also use these variables to re-express the fields in 5^^*. We will then have 


^(3) + /s:(3) = ^ J ^3^ ix,t)Ux,t) + X ^ ^ 

a fiyX 

In Appendix B, we show that 






We define 




(Or COc 

p = I ; sin a 


16;r nokT 


|9£/./3<oj 
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and obtain 


AT<^) = 2^ ^IJsJlJt cos <t> 

where 


<p ss 01^^ Bt — Os 
We will also define 

ro^P 4Jl 

The factor of 2 comes from the fact that the growth rate for the action is twice the growth rate for 
the fields, as computed by Shukla et al. (1983). 

Going back to the action + and varying with respect to we get 

-^ = ^ + V, • ? /<, = -2^ A sin 4) (7.1.1) 

We also get: 



From 7.1.1 we can see that — = 0 unless (9^ + ^v + ^a.) not vanish when the bar indicates a 

at 

lime average. This gives us the resonance conditions 


(Ofy- and = ^ (7.1.1a) 

If one wave, initially has large amplitude and the other two 7^, have very small ampli- 
tude, then initially we may take as almost constant. Furthermore, the relative phase <P will 
quickly evolve so that sind>= 1, and so we use sin<I>= 1 in what follows. (For this and a discus- 
sion of the general time development of such a system of equations see Meiss, 1979, and Berg- 
mann, 1985.) 

We define = h/tJ a-y = obtain 



+ V, 


S\ 


doi 

17 


= ro^2 


3 ^ 3^2 


= ro(i\ 
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If waves 1 and 2 are damped, we can put in the damping heuristically to obtain the WKB equa- 






which were first obtained by Rosenbluth (1972). 


(7.1.2a) 


{lA,2b) 


In the cases where the damping is of kinetic origin (ion Landau damping for example), this 
would have come out of the solution of e{k, ct)) = 0. In the case of the transverse waves, however, 
where the damping is dominated by collisional effects, it is not deriveable from an action principle. 

A brief word about the damping rates of the ion-acoustic and electromagnetic waves is in 
order. For the transverse waves we use the Dawson-Oberman impedance 

r'r.2lnA-l . . . .3 

< — ^ — where A = 4;wA^ 

<»>pe Vl8n:A 

{Dawson. 1968). This may actually be as much as a factor of 2 larger due to enhancement from 
ion density perturbations in the solar wind. 


— — ~ 2.5 X for the 11 March 1979 event. 


(7.1.3) 


For the ion-acoustic wave, the damping rate as determined by Lin et al. (1986), is 


rs = A(o 


(7.1.4) 


due to Landau damping of the ion-acoustic waves. For the 11 March 79 event T^/T, - 5. 

The damping of ion-acoustic waves leads to a dissipative “localization” in space. The dissi- 
pation length will be 


g _ ^5 ^ 1 

Ys -1 
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7.2. Space-Time Evolution 


If we assume that a^{xj)=-a ^ e- and a2(A:,/)= 2 t 2 and substitute these into 

lA.lb, then {co- K'Vx^y\)'^\-yo^2^ («^-^*V2 + r2)^2 = To 2 ri, which can be combined to 

yield the “dispcrson relation” 


(o) - K'Vi + coi)(co-^- V2 + 72)” 75 = 0 • (7.2,1) 

The spatial and temporal development of a 3-wave system governed by 7.2.1 have been 
thoroughly studied and reviewed by A. Bers (1983) whose analysis we follow in the remainder of 
this section. In order to simplify the analysis, we first transform to a frame moving with velocity 


V such that Vj = vJ+l/, V2-V2 + V, and a>' s co-^ K where V is chosen such that V|,V 2 are 
(anti-) parallel. We will take wave 1 to be the ion-acoustic daughter wave and wave 2 to be the 
transverse electromagnetic daughter wave. The general vector relationship of and are shown 
in 7.2.2, along with our choice for V, which will reduce the problem to one-dimension with the 
vector relation shown in 7.2.3. 


Si 


7.2.2 


V ' 


7.2.3 




In the frame of reference shown in 7.2.3, one obtains the “dispersion relation” 


(co' - K Yj + 7 ,) (co' - Kv's + yi)-^=0 

If 7| = 72 = 0, this would be an absolute instability (positive growth rate for a fixed position in 
space) with the growth rate 


2 1 vj Vs 1 

^ ^ 7 - 7o . 

I Vj I + I Vs ^ 

There are, however, two thresholds associated with and 72 ^>^0: 


1) 75 >7i 72 for there to be an instability at all, and 


(7.2.4) 


72 


^ [yi I v\ I + I vj- 1 ]• 

2 ) y~> 7 -; for the instability to be absolute. 


(7.2,5) 


4 I V 5 VY I 

If (7.2.4) is satisfied, but (7.2.5) is not, then the instability is convective, meaning that at a fixed 
point, X, in space a (jc ,0 — > 0 as f but there is a moving pulse that will grow in time. 

The growth rate of a convective pulse moving with velocity u parallel to Vs, and 7.2.3 is 
given by 


„ 2 [(Vr-«)(Vj + M)]'^ yi(M + V5) + y2(V7--M) 

^ yo — 


(7.2.6) 


VY + V5 V5 + Vx 

For the March 11 case 70 = 1.8/s, as was determined by Lin et aL (1986). Using 7.1.3 and 
7.1.4, 75 >ys yr so the process is unstable. On the other hand it is the case in this problem that 
Yt '’i fs 'Y- and so. 


{■/ I V5 I + 75 I v't I f 


— =4x10® 

4v; 


4 I Vs lY I 

Therefore the condition (7.2.5) is not satisfied by 6 orders of magnitude. When yj — > 0, (7.2.6) 
becomes 


2 [(v^-iO (Vs + «)]''- yi(vr-«) 

"0 = ^ ^ yo : — 7“ 

VY V 5 v ’5 "h Vj* 

It is easy to show that the fastest growing pulse will be when 


(vr~u) = 


vry? 

7s 


which gives 


o)„ = 


7s 


The growth length for the 1 1 March 79 event will be 


Ls — - = = 7 X 10 ^ km . 

^ (^0 yo 
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This presents a distinct problem for growth. An ideal, non-turbulent solar wind plasma 
would have a density gradient (due to conservation of the solar wind flux) 


2n 


$Q) 


where R is the distance from the sun to the point of interest (1 AU). The eikonal equation for the 
wave vector h of a wave propagating in media with inhomogeneities with a scale length much 
longer than the wavelength of the wave is: 


dri _ JL Y ^ /} {h ♦ Vv) 
dl 


0^2J) 


V 

where v is the group velocity of the wave {Landau and Lifshitz, voL 6, para. 66). Since ki is 
approximately parallel to the density gradient in the solar wind, vj will be transverse to it. The 

dominant term in 7.2.7 is therefore — Vv, and we find that 

V 


4 = 1.6 xlO“*/km . 


dl c''-k- R 


In one growth length the wave will be displaced a distance along the density gradient given by 


-4 = 3920 km. 

dl 2 

The resonance condition for growth is such that growth will cease when v^AA' =yo. But 
V- Ait = \> Vfc A/? = -^ AR=2 

* * dl (Op 

and so the transverse wave will be refracted out of the resonance region within oen growth length. 
The consequence of this is that because of the low level of the Langmuir pump and the high damp- 
ing rate of the ion-acoustic daughter wave, the instability is convective instead of absolute, but 
because of refractive effects the convective instability is suppressed. 

By contrast, we apply the same considerations to the Whelan and Stenzel (1985) experiment. 
Using their results we estimate 


7o = 3x 10^/sec^, = 1 X 10‘Vsec“ 

vj = 4.5 X 10’ cm/sec and 
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\’5 = 2.6 X 10^ cm/sec 


We see that To » yr 


vr 

4vc 


and therefore in this case the absolute threshold is greatly exceeded. 


If the discussion in this section has so far cast doube that a 3-wave process is involved, the 
following argument gives strong support for it. If we write down eq. 7.1.1 for each wave in the 
type III process, assuming no wave damping, we get 


— ^ = sin CD 

at 


dt 

dt 


— 2/3 sin<I> 


dJ: 


- = 2/3 VAA/rSind> 


dt 


Adding the first two gives us 


dt 


(Jl^Js) = 0 


(7.2.8) 


This is one of the Manley-Rowe relations (Sagdeev and Galeev, 1969). Again ignoring damping 
and convective effects, we would expect from eq. 7.2.8 that at saturation of the 3-wave decay 
=75. Using the standard expression for £5 (/:, <o) and Ci (k, o>) gives us 


1 




Jl 




= 1.02 


when we substitute parameters for the 11 March 1979 event. Bergman (1985) did an extensive 
numerical study of the time evolution 3-wave systems when one of the daughters is damped and 
found that the damped daughter wave and pump wave will oscillate in amplitude 180® out of phase 
with each other, between zero and nearly equal maximum amplitude, which, however, decay with 
the damping rate of the damped daughter wave. Actually, this argument ignores another effect 
which is that the electron beam will replenish the Langmuir pump wave as the latter is depleted by 
the 3-wave decay. 


Because all other factors seem to support the 3-wave process we next examine how density 
turbulence modifies the convective process and what effect this modification would have. We begin 
with a review of the considerable effort that has been made to understand the effect of spatial 



75 


plasma inhomogeneities on parametric instabilities, much of which has been directed towards 
understanding laser-induced plasma instabilities for laser fusion applications. 

The usual approach to this problem begins with the following equations for the slowly evolv- 
ing amplitudes 


da\ 

( ^ 



( \ 


, dx 

= yo^2 

exp 


dt 

'da;' 
dx J 

= ro^i 

exp 1 



where K - {Rosenbluth, 1972; Nishikawa, 1968). For a linear density profile such that 

t 

K{x)- K' (a*-ao), Rosenbliith concluded that for V] V2<0, when both waves propagate along the 
density gradient, the gradient leads to convective growth of the instability. DuBois et al. (1974) 
showed, however, that this effect was due to the unphysical assumption of an infinite length plasma 
with constant density gradient, and finite lengths lead to temporal growth for normal modes of the 
system. 

Klein et al. (1973) considered the case of the Raman side-scatter instability, the geometry of 
which is similar to the type III problem. A laser pump propagating along the density gradient 
decays into a Langmuir wave propagating approximately along the density gradient and another 
electromagnetic wave propagating transverse to the density gradient. The authors did a 2-D numer- 
ical simulation and discovered that this instability is not stabilized by a strong density gradient as is 
the Raman backscatter instability. This was attributed to the side-scattered wave remaining in the 
resonance region until refracted out, allowing for a much longer growth length. 

Further theoretical progress on the Raman backscatter instability was made by Drake et al. 
(1973), after Forsliind et al. (1973) discovered through computer simulation that a temporally 
growing mode (absolute instability) existed in the Raman backscatter problem. Drake et al. 
pointed out that the beating of the pump and the backscattered wave caused a bunching of electrons 
which would partially reflect the backscattered wave back into the region of the instability. Tem- 
poral growth became possible for waves which were quantized by the well formed between the two 
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turning points. 

M. Mostrom (1975) considered this effect for the Raman side-scatter problem using WKB 
techniques. The conclusion was that the reflection due to the inhomogeneities produced by the 
beating of the pump and side-scattered wave could balance the effect of refraction, thereby allowing 
for absolute rather than convective growth. 

The first study of the effect of random inhomogeneities on parametric instabilities was done 
by Tamoikin and Fainshtein (1972). The effect of random inhomogeneities considered in this study 
is to cause the phases of the three waves to fluctuate randomly, and to be scattered. The major 
effect is shown to be an exponential damping of the pump wave, due to scattering by the inhomo- 
geneities. This depletion of the pump damps the parametric decay. 

Nicholson and Kaufman (1974) studied the effect of random density fluctuations superim- 
posed on a linear density profile. The suppression of a temporally growing mode, as discussed by 
Rosenbluth (1972) and mentioned above, was due to phase cancellations which, as Nicholson and 
Kaufman showed, is destroyed by random fluctuations. They therefore showed that a temporally 
growing mode can occur. Nicholson (1976) also studied the effect of a sinusoidal density modula- 
tion on parametric decays. The effect of this density modulation was discovered to suppress 
growth rates, but a large density modulation was required to provide substantial suppression. 

The effect of random wave phase caused by time varying random density fluctuations was 
studied by Laval et al. (1976) using the Bouret approximation, who concluded that the effect was to 
suppress growth. 

In summary, the major effects that have been considered are the effects of inhomogeneity in 
the phase relationship between the waves. Only in the studies of Drake and Mostrom has the effect 
of reflection and partial trapping by inhomogeneities been considered, yet for the cases studied the 
daughter waves propagate near the critical density where even small inhomogeneities will have 
major effects. 

For long wavelength large amplitude density fluctuations, on the other hand, waves near criti- 
cal density will be completely trapped. Although such a trapped region will be very irregular in 
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shape, leading to chaotic ray trajectories, a large body of work has been done which shows that in 
such cases, waves will in fact be quantized. This was studied by McDonald (1983). There is also 
a large body of work in quantum chemistry studying the quantization of classically chaotic systems 
(see Casati, 1983, and references therein). 

The effect of wave-trapping on a parametric instability is two-fold. The first is that since 
there is no wave convection, the instability will be temporally growing rather than spatially grow- 
ing. The second is that the trapped wave will be quantized. The discrete spectrum for k and co 
will limit the possible waves in the instability and may slow or inhibit the growth due to phase 
mismatches. The effect of smaller amplitude shorter scale length fluctuations is not so easily 
treated in this scheme, and for this situation we return to localization effects. 


7.3. Localization and Three- Wave Processes 

We are nearly at the end of the story. The effect of wave localization on the 3-wave process 
is to prevent wave convection out of the region of interaction non-dissipatively. Localization of the 
Langmuir waves is unnecessary for this, but as we saw in section 7.2, prevention of convection of 
Che transverse wave is necessary. In order to make theoretical progress we assume, again, that once 
localization effects on wave propagation are included, we can use averaged plasma properties to 
compute coupling constants. We assume therefore that the Langmuir and transverse waves are 
exponentially localized so that we may use a modified eikonal representation of the waves as fol- 
lows: 


Er = Eje (k,x) e‘^ + cc 

= e'^‘- + cc 

Then, 


5<-' + 5<^» = jd^^dt [eL(kL,x) E[ + Dt (kT,x) Ef 

+ es(ks,x)Ei(x,t) 

+ lijd^xdt sina(.v)£f. £^ J 


Xt\ 


+ CC 
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where 




npe' 

m^kTM o>l a>5 


as determined in Appendix B. Again 
h '^T 


sin a(x) a 


IJfcLl 


where in this case it may be a function of position. Take variations 




dOj * dt dco 


^ V L -r 


_i_ . 0 = -^ +i^ smaix)Ei. 

dSj dt dco 


^ V L -r 


— :0 = -^Ei-ip sind(x)ELETEs e e 

60s dt dco 


If we define the total wave action 




jd^x \E[\ 

[acoj 


r^l _ de 

.3 dco 


2 


J,.U>. I ^ , i£iL ^ 

^ J d<o d(o 

kt 


where we have assumed that ‘ , and therefore that the localization i 

Then 


and 


-VS'- 


Kiix-Xi^l 




{dDIdco) 


T -~Kr\X — Xri 

— e 


(7.3.1a) 


(7.3.117) 


(7.3.1c) 


isotropic. 
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Since the ion-acoustic wave convects and is damped, we must retain the action density 
If we integrate 7.3.1a-c 


-^J - 




Ir -3 fc .3 
^Ti 


— Iff Ijc - Xj - 1 


X - Yl Jl 




'Sl (Be,fd<o)(iesZ)(dDr/d<o) sin«(a:)cos e 


^3 „3 


— I X — X£, I t — JCf I 


X aJJi JfJ^^x) 


~YtJt 


ijsU) = 


2p 


■4 


Kl Kj 


{dEJd(o) (dE^/^o)) QDjldco) 


sina(x)cos(t>(x) e 


— if£, I X — x£ I ^”*‘T J Jf — JTr I 


X ^Jji^ ‘/r/sC'^) 


-ysA(-v) 


As before, we consider the case where at / = 0, is large and £ 7 -, £5 are at fluctuation levels, so 
that initially we can take £ 7 , a constant. Then we use the definition 






(jdejhoi) (des/Boi) (dDf/^oi) 

With this definition yo will be of the same form as in section 7.1 where £7. is now the peak value 
of the electric field. Again we have 

■^ + )v| ^Jj^ = Yo‘>Jk^ jd^x |sina(.v)cos^(.v)e"‘"‘-'^"^'-' ■ 


^ d ' 


V Js ) = 7o sin a(.v ) cos 0 (.r ) e ^ 


and 
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Combining these gives 




( \ 



+ Js 


sin a(jc)cos <p (jc) 


^ \ 

- jd^x' sina(;c')cos^(jcO^”’^^*^*’''^' . 

Although we could attempt to integrate this numerically over assumed forms for the localized 
waves, we can make some estimates as follows. 

Experimentally, it appears that the envelope of the ion-acoustic wave is similar to that of the 
Langmuir wave. Furthermore we assume 


sin a(A') - 1 and cos <^(a^) ~ 1 . 
Therefore, we assume we can write 

We make the WKB assumption 

Then 


{(o + Yt) (oi + y^) = kI ^jd'x 

Let us further assume that Kf and that I - Xj I 


«K— . Then 

Kt 


(co + 7r)(^ + rs) = 




To 


which gives us the requirement for positive growth that 


To > 


2kt 


1 3 


YtYs 


For the 1 1 March event 


(7.3.2) 
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70 km 




1 


500 km 


— -2.5x10-"; 

Ys_ 


0>5 
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The instability requirement (7.3.2) is >9 x 10“Vsec^. The actual value of 

75 =(1.8/sec)“ = 3.24/sec" and we see that the instability threshold is exceeded by more than 3 ord- 
ers of magnitude. The effects of phase decorrelations, etc., might be estimated at 1 to 2 orders of 
magnitude. 

The following considerations are of possible significance: 

i) Because the Langmuir localization length is much less than a transverse wavelength, wave 
number matching is not significant. We should merely have ks = k^^ On the other hand, fre- 
quency matching is significant. It is necessary that the frequency mismatch be less than the 
growth rate Yq. The frequency matching will be determined by the frequency levels of local- 
ized states. 

ii) Other than assuming that k^ is perpendicular to ki, we have not really considered the effects 
of the 3-dimensionality of the problem. Although a one-dimensional model may be adequate 
for the Langmuir waves, the transverse wave is certainly a 3-dimensional problem, and the 
localization length predictions from the 1-dimensional theory are certainly shorter than in real- 
ity, Furthermore, the electromagnetic wave is a vector wave and consideration of the effects 
of scattering on the wave polarization is important. 
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8« Concluding Remarks 

We have argued that wave localization effects can be important whenever a wave propagates 
near the wave cut-off in a plasma with superimposed density turbulence. The effect of localization 
will be to stop propagation and limit wave spreading to a velocity approximately that of the density 
fluctuations. Furthermore, the confinement of the wave and possible daughter waves can change 
the nature and threshold of plasma instabilities. 

We have discussed beam-plasma interactions at length because this problem is easily treated 
by one-dimensional simulations. The localization of Langmuir waves in itself would seem to be a 
curiosity with little importance other than proving that density fluctuations do not suppress the 
beam-plasma instability. On the other hand the localization of the transverse waves allows an ins- 
tability to reach threshold which otherwise would not occur. 

We have only simulated those parts of the problem amenable to a one-dimensional treatment. 
To treat the full problem, including 3-wave decay and non-linear effect in even two dimensions, by 
computer simulation would be very difficult. Because effects happen on both electron and ion time 
scales, a two component model of the background fluid would be needed, this leads unfortunately 
to charge separation and noise for warm fluids. A large number of grid points would be needed 
since the system needs to be large enough to include several electromagnetic wave wavelengths. 
When beam particles are included one obtains a memory requirement of -4x10^ words. This 
would be prohibitively expensive in Cray time. 

Although the existence of wave localization in 3-dimensions is well established theoretically 
for fixed scatterers, the problem of mobile scaiterers in 3-dimensions is not settled. This problem 
is amenable to study using an extension of the methods used in Chapter 5 on, say, a 
lOOx 100 X 100 lattice. In this case, only waves localized to a small number of lattice points could 
be studied. 

Within the context of the type III problem, the arguments we present are consistent with 
observations at 1 AU. We make no claim as to their importance or validity in other regions, in the 


83 


solar corona, for example, where strong turbulence effects may in fact be important. On the other 
hand it seems improbable that if the solar wind is turbulent at 1 AU that is is not also turbulent in 
other and perhaps all regions closer to the sun. In this case, localization will still occur for small 
amplitude waves. For waves of large amplitude, strong turbulence effects including soliton forma- 
tion and collapse are probably of more importance, although wave localization may cause a faster 
onset of these processes. 

A final word for those who would argue that wave localization is unnecessary for understand- 
ing the behavior of the Langmuir waves. To a certain extent, we would agree. As the scale length, 
L, of the shortest wavelength density fluctuations increases to the point that L semi-classical 
or WKB techniques become valid. As McDonald (1983) showed, the closed, quantized trajectories 
are of special significance. Since we are ultimately concerned with wave behavior, however, we 
argue that one cannot just use ray trajectories and ignore quantization effects. This is the import of 
the work by Bezzerides (1986), discussed in the introduction. Even if, because of particular aniso- 
tropies in the turbulence, most trajectories are not trapped, some will be, and when closed and 
quantized, lead to long-lived quasi-stationary states. Once the Langmuir waves are confined and 
quantized, whether through localization or semi-classical effects, much of the discussion of Chapter 
6 holds. Because of the length scales involved, the properties of the transverse waves can only be 
understood through the properties of the full wave equation in random media, i.e., localization 
theory, since dispersive effects are as important as refractive effects. The advantage of using local- 
ization theory, on the other hand, is that one can be very ignorant about the turbulence involved 
and still make very powerful statements about the consequences. 


84 


Appendix A. Simulation Techniques and Code Listing 

The simulation codes developed for this study are based on standard plasma simulation tech- 
niques as discussed in Birdsall and Langdon (1985), which is the reference for the numerical tech- 
niques discussed below, and some of the code is taken directly from the program ESI by Bruce 
Langdon. 

For the localization studies we have used two programs, BEAMLOC, and PERLOC. BEAM- 
LOC is used for the beam-plasma interaction, and is listed following this discussion. PERLOC was 
used for the Gaussian wave packet propagation, and is similar to BEAMLOC except that the 
unnecessary particle mover and field solver were eliminated, resulting in a very much faster run- 
ning program. 

In the remainder of this section we discuss BEAMLOC; however, the discussion of PERLOC, 
suitably restricted, would be the same. 

Our techniques for simulating the background plasma is to use a linearized one-fiuid plasma. 
Our rationale for this is the following. 

In this problem we wish to simulate the effect of a very weak beam with physical density 
10“^. If we tried to use particles, in order to observe the fields that arrive from the beam 

plasma 

plasma interaction, we would certainly expect that the statistical fluctuations in the background den- 
sity would be significantly less than the beam density. 

Let = number of background simulation particles in the box A.r; n^= number of beam 
simulation particles. The charge of a simulation particle is 


4 k (e/m) Up 

ie/m) is established by the particle type. 

If we require that the fluctuations in the background particles produce fields significantly less 
than those arising from the beam-plasma interaction, then 
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£ = fluctuating electric field 
n = fluctuating density 



1 ^ = -^ 1 

k elm k 


V E B Anen 



(^/w) 


•y 



We require 


We get 

* "fc 1 

— ^ Ax <«: 

{€im)n {elm) k 


rip (^'A.v)- 



2 




It is quite impossible to use particles for the background in this case. 

A further problem, which also restricts the ability to use a two-component fluid for the back- 
ground is that since we have significant density gradients because of our fluctuating background 
fields, we will get strong electric fields due to charge separation unless we have Ax<i:Xp. How- 
ever, for our problems we want to follow several hundred wavelengths and want 100 Debye lengths 
per wavelength. In order to satisfy all of this, we should then require >10^ grid points. 

By using the equation of motion 


^ = -a>;(x)v^ + 3vr ^ 

Or ^ Or 

for the turbulent background plasma we gain in simulation quietness and isolation of the effect we 
are interested in (the effect of density fluctuation on wave propagation). On the other hand, we 
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completely eliminate any effects which are non-linear in the background plasma. This could be 
significant since the density fluctuations are large. 

SPERLOC uses the following equations. We assume a grid with + 1 grid points, equally 
spaced over a length L giving A.r = — . Quantities defined on the grid are advanced by time steps 
A/. 

The equations of motion we solve are: 

Particles: 


d'-.x, ^ 

dr >n 

Wbeam (.Cl /) = X ^ 


Fluid: 


Define v'(v,<) = -3 [«a(-v,/) - 


V.3v, 


where co;(.r)s 


Anui (x)e^ 




Field: 


V • £ = Anp where p = + ¥ 

Initial conditions: 


V^U,0) = 0 


co;U) = 


1± 


"0 


where 5ri{x) is a randomly distributed variable 


x,(t =0) is uniformly distributed in space 
\r{f =0) = drifting maxwellian distribution 


87 


We use a time centered leap-frog scheme to solve these equations. We define the quantities: 

= AL y (t + ^/ 2 )At) the velocity for the y th particle 

^ U = ;A^; t = (n + '/.)A/) 

x"sx,(t = nAt) the position of the yth particle 
xf/j s i//(,r = jAx; t = nAt) 
oiji = oil(x = iAx) 

We first solve for the electric fields in Fourier space from the following. 


P.— Vf + P, 'beam 


A-£," = 4;rp; 

and then advance quantities in time: 


,, —ejAt 

J T?n . ..n 

V/ V,- 


rt-Vi 


\l/r^ = XI/r + At 
.v,"+' = .r;+vf+'^A/ 

We then find p!"*', then E”*' and so on. The advantages of a time centered leap-frog scheme to 
solve differential equations which are 2nd order in time is that the solutions are of accuracy 0 (At-). 
Plots: 

In addition to plotting the quantities x(f" and we are also interested in the “energy den- 
sity,” which we define as 

r-n 1 •> , 1 • . 3 , [V/+I - ¥i-\]' 

E. - J V? t 2 V, • 2 V,- 

This is the Hamiltonian density which will generate the equations of motion for xy and xj/. Plotting 
this has the advantage that the quantity is phase independent, unlike x(f and xj/. 
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Model for background density fluctuations 

In most cases the model we use for density fluctuations is that of Escande and Souillard, 

which is that co^{x) is a step function, constant on a specified number of gridpoints, which is 4, 
8, or 16 in most of the simulations presented. The value at each step is generated by a random 
number generator and is uniform on the interval. 


We show the power spectrum of density fluctuations in this model in Figure B.2. In this 
figure / = 314, n^ = 512. According to this, 

dx 8n{x) = ^ f * 

^ Jo jA 

This gives 



M} 

' sin^ kljl ' 

»'=0 




For moderate k, the sum vanishes unless j = j' and so 

sin“ kljl 
ikljlf 

As ^ 0 however. 


j 


I Sn^ I “ — > 






= 0 


The peak occurs for kl^ = .55 


512 


X lOO/r 


= 1.227 and decreases approximately as ^ ^ up 


to kl^ - Ik, The detailed structure, however, is of course a specific manifestation of the statistics of 
the ensemble, since the sample is finite, and it is this that gives rise to the localization. 

Notes on units: 

We use the following convention for units in this program. We choose units of time so that 
the background plasma frequency co^^^=l. We choose units of mass so that e/m = l. Once we 
specify the number of simulation particles /?bcam ^ t)Ox of length Ar, we have set the charge of a 


89 


beam particle 

2 

_ beam 

4;r(e//w ) «be„„ ' 

Typically we take 

2 

_ ^bcam _ g ^ |Q-6 
< " " 

As it is, we find that 4,096 grid points is the largest number necessary, and 8,192 particles 
gives us two simulation particles per grid point. 

V, 

A typical length of the system is L -314, and v, ~ .05, With co = 1, Xp = — = .05. With 

4,096 the number of gridpoints, we have Axlkp=^ 1.5. As none of the physics of our simula- 
tion relies on being able to resolve a Debye length this is sufficient. Also a typical beam velocity 
will be V/, = 1.0. Then 

A'o= — = 1.0 A, -4^ =6.28 — = 82, 

V* ^0 Ax 

therefore a typical wavelength is very well resolved. 

The typical CRAY run time for BEAMLOC with /j^= 4,096, 8,192, /j^= 4,000 was 10 

minutes. The typical run time for PERLOC with =4,096, rtf = 80,000 was 4 minutes. 
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Appendix B. Computation of the 3-Wave Coupling Constant 

In this appendix we calculate the 3-wave coupling constant between a Langmuir wave, an 
ion-acoustic wave and a transverse electromagnetic wave, used in Chapter 7. We begin with 
eq. 7.1.1, 

(B.l) 

In this particular case we note that will contain no term of the form since 

- Ml - and the only wave with non-vanishing A in our choice of gauge is the transverse wave. 
Let us write 




i m,c 

where a is a wave index and = [//(')]’ (complex conjugate), and where we have defined the 
“non-relativistic” 4-vectors 


A = (0, Aj) for the transverse wave 
i ~ 

A = (<?a,0) for the ion-acoustic or Langmuir waves 


and pi = imiC,m,Vi). 

Then = jdt 

where = -cu„. Eq. B.l becomes 

I 

i 

iS:<3)=l£{ W,{VV, //<*>) ) 
a 

U* we use the Vlasov approximation 

/ (.v,p ) = (.V - .V,) 5(p- Pi) 

i 

then 

dt J (vt-,, } = X jd^.xd^p ^5 (.V -a:,) 5 (/>-/>,) {n>, (w, //<’>) ) 

= -jd^xdt d^prtogip) j i w, ) 

I 
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Integration by parts gives 


= T Z ^ V "0 { w , g (p ) } { w . //i*> ) 


= \'^ld\xd^p no 


dw dg 
[ dp j 


{M',//i'>) 


If g(/7) is a Maxwellian distribution, then 


dp mlcT 
and we find 


g(P) 


S'»= ]■ S j(;,) I 


We take 


w= > — 

a 


which gives 


9w V* // \ 

'’ SJ-I-A-**-'’* 




= I =I 

Combining the above gives 


i9„ i9„ 


K^^^^\\d\xd^p 


nogip) 

^ \a xa p 

3 J ^ mkT 

Now, we note that for :» Vj- , 


a 


z 

/Jy 






i9a 


V = (Hl 


V* 


= 0)^ 


0)7* - Xy • V = <0£ since <«: it/. 

and we assume that we may use <05 - it^ v — • v when integrated over the thermal distribution. 


We write to represent the sum over all permutations of ajiy where one of the waves is 


a Langmuir wave, one ion>acoustic, and one transverse. 
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ks ■ p ki • £j ^ p • (.ki • kj) ^ (kf ■ p){ki • k^) 1 ^ 1 

Os Ol n9? ®T flf er 


ki P ks ■ St ^ (/7 • er) ks • kj ^ ki p k$ -ej ^ kj-ks 
Oi mdr mOf 


ks -p ki -B t ^ p -Cr (^t, ‘^d 

0c &T ^ flj 


Since ki = ks and 9[^ = 9t'»8s^ the underlined terms are approximately equal and are dominant. 


cancelling the —. If we define 


ks'h 


5<^) = jdr d^x d^p sin a ** *^^ + c.c. 

J kT me 


r X ^0 k$ . - — »—• i(o,-0,~Or) 

^{dtd^x— — sin a At- e 

J Jtr me a>T, 


We evaluate the quantity /3 defined as 


'*0 ( 1 


as follows: 


kT 1 me 1 o>i 


Noting that 


= . "r 7 “i- 7 

Ef — i AfBi Af 

c c 

Es = iksPs 


El = iki <t>i 

and using the definitions 


4/r dct> coi 
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/s = 


^ ^ " de 

4tc 3o> 


= l£c 


2(o} 


cos 


where we have used the standard expression for the dielectric response functions and Then 
we find 





16;r no kT 


sin a 


The factor 


16k riokT 

is the same as that derived by Shukla and Yu (1983) using fluid equations. We can also write p in 


P^lPs^t ~ 


m 


iripe 

^kTM a>l CO5 


ElEs Ej 


the form 
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Figure 2.1 Electron flux and 
electric field measurements 
(from Lin et. al., 1986) 
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Figure 2>2 Reduced parallel electron 
distribution function, (from Lin et. al,, 
1986) 
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